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KXXIX. BAND ERSTES HEFT , 1960 


Kinige Betrachtungen zu den Grundlagen der Umformtechnik * 


Von Th. Lehmann 


1. Einleitung. Die Umformtechnik umfaB8t alle Fertigungsverfahren!, bei denen durch bleibende 
Forminderungen im festen Aggregatzustand aus einer Ausgangsform eine neue Form, die Endform, 
entsteht. Der Ablauf des Umformvorganges ist durch die zeitliche Folge der Augenblicksformen 
gekennzeichnet. UmfaBt ein Umformvorgang mehrere Stufen, so erhalt man Zwischenformen als 
Ergebnisse der einzelnen Umformstufen. 


_ Die Umformung ist ein irreversibler Vorgang, und zwar nicht nur im thermodynamischen 
Sinne ; denn jede Umformung hinterla8t im Werkstiick Veranderungen, die praktisch nicht mehr 
oder zumindest nicht mehr vollstandig riickgingig zu machen sind. Hier sollen aber nur solche 


Erscheinungen betrachtet werden, die phanomenologisch im Sinne der Physik der Kontinua er- 


faBbar sind. Kinzelheiten der Vorgange innerhalb der Kristallite, also Gefiigeanderungen, Ver- 
anderungen im Kristallgitter usw. oder, anders gesagt, alle Vorgange, die in unmittelbarem Zu- 
sammenhang mit der atomaren, molekularen bzw. kristallinen Struktur des Stoffes stehen, scheiden 
aus dem Kreis der Betrachtungen aus. 


Phanomenologisch gesehen handelt es sich bei Umformvorgingen um ,,Zustandsdnderungen“ 
des Werkstoffes. Zustandsanderungen, d.h. Anderungen einer ,,Zustandsfunktion (abhangige 
Variable), ergeben sich aus Anderungen der ,,Zustandsgréfen“ (unabhangige Variable). Zustands- 
nktion und Zustandsgré8en sind durch ,,Zustandsgleichungen“ miteinander verkniipft. Die 
Jefinition von ZustandsgréfBen ist an die Forderung gebunden, da8 die Zustandsfunktionen ein- 
Jeutig von den ZustandsgréBen abhingen miissen, was in der sogenannten Integrabilitatsbe- 
ingung der Zustandsfunktionen zum Ausdruck kommt? %. Deshalb kénnen die Variablen eines 
organgs nicht in beliebiger Weise als ZustandsgréBen eingefihrt werden. Hine weitere Forderung 
= daB ZustandsgréBen-stets unabhangig vom Koordinatensystem, d.h. invariante Groen sein 
lissen. 


Die Frage nach der Definition von ZustandsgréBen bei Umformvorgangen ist bisher nicht 
ligemeingiiltig geklart; sie ist wohl auch noch nicht in so allgemeiner Form gestellt worden, wenn 
an auch den Zusammenhangen mit der Thermodynamik schon einige Betrachtungen gewidmet 
t4,5, Es ware sicher niitzlich, dieser Frage nachzugehen; denn in engem Zusammenhang damit 
ehen praktisch wichtige Fragen, wie z. B. die Frage nach der ertragbaren Beanspruchung, die ja 
eine Funktion der ZustandsgréBen sein kann. In dieser Arbeit kann dieser ganze Fragenkreis 
och nur kurz gestreift werden. 


enden Formanderungen zu tun. AuSerdem hat man im allgemeinen damit zu rechnen, daB 
ie Werkstoffeigenschaften mit der Formanderung verandern. Diesen Gegebenheiten wird man 
ten dadurch gerecht, da8 man zur Beschreibung der Vorginge ein kérperfestes Koordinaten- 
em verwendet. Ausnahmen bilden gewisse stationdre oder quasistationare Vorginge, die sich 


MG 


Die nachstehenden Ausfiihrungen sind der von der Fakultat fiir Maschinenwesen der Technischen Hoch- 
Hannover genehmigten Habilitationsschrift des Verfassers (Berichter: Professor Dr.-Ing. E. Pestel und: 
ofessor Dr.-Ing. O. Kienzle) entnommen. Mit Riicksicht auf die darin ebenfalls enthaltenen Betrachtungen 
rfassung und Ordnung technischer Umformvorgange erhielt die gesamte Arbeit den vorstehenden Titel. 
Fertigungsverfahren sind dadurch gekennzeichnet, da ihre Ergebnisse stets feste K6rper mit definierten 
sind; Zubereitungsverfahren beziehen sich dagegen stets auf ,,formlose™ Korper, d. h. auf solche ohne 
rte Formen, wie z. B. Gase, Fliissigkeiten, Pulver, Granulate, formlos erstarrte Flissigkeiten usw. 
R. Becker, Theorie der Warme, Berlin/Gottingen/Heidelberg 1955. 
Hausen, Brennstoff-Warme-—Kraft 2 (1950) S. 1, 43, 77. 

VU. Freudenthal, Inelastisches Verhalten von Werkstoffen, Berlin 1955. 
Ziegler, Ing.-Arch. 25 (1957) S. 58. 


ei der Betrachtung bleibender Formanderungen hat man es in der Regel nicht mehr mit klein- . 
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auch in ortsfesten Koordinatensystemen und mit Methoden verfolgen lassen, wie sie in der Str6- 
mungstechnik und in der klassischen Plastizitatstheorie 12 3 4 5 6 78 iiblich sind. a) 

Bei den folgenden Untersuchungen wird die analytische Tensorschreibweise nach DIN 1303 | 
verwendet, unter Einschlu8 der Einsteinschen Summationskonvention. Im iibrigen lehnt sich i) 
die Art der Darstellung an A. E. Green und W. Zerna™ an und baut wie diese auf dem Ricci-Kalkil i 


auf. 


2. Koordinatensystem und MaBtensor. Das kérperfeste Koordinatensystem sei mit x' bezeichnet. | 
Von ihm wird neben gewissen Stetigkeitsforderungen zur Sicherung hinreichender Differentiier- | 
barkeit nur vorausgesetzt, daB jedes Koordinatentripel einen Punkt des Korpers umkehrbar ein- | 
deutig festlegt. Im iibrigen ist das Koordinatensystem beliebig. Es macht alle Formanderungen | 
des Kérpers mit, so da8 jeder Kérperpunkt seine Koordinaten bei allen Formanderungen und Be- 
wegungen beibehalt. Bildet man das kérperfeste Koordinatensystem in einem ortsfesten, carte- | 
sischen Koordinatensystem (y') ab, so hat man im allgemeinen einen nichtlinearen Zusammenhang | 
zwischen den Koordinaten, der auch noch von der Zeit t abhangt: 


yi =y' (x9), ai = ai (y*, 2). (1) | 
Die Koordinatendifferentiale transformieren sich dagegen linear: | 
dy =e ae; da = ciady". (2) | 

Die Transformationsmatrizen 
: ayi pepe | 
bee Tao? Cre dyk (3) | 


sind im allgemeinen Funktionen des Ortes und der Zeit. 
Die Metrik des kérperfesten Koordinatensystems ist durch den symmetrischen MaBtensor fest- 
elegt!8; 
6 8 es 5 heed 1 oe eal te Soe) (4) } 
Sie SF os Be Oe. Ce 


GréBen, die sich auf den unverformten Zustand des Kérpers beziehen, werden durch eine tiber- 
gesetzte © gekennzeichnet, z. B. ik oder ch. 


3. Formanderungstensor. a) Definition des Formanderungstensors. Aus der Beziehung | 


Bik = Bim 8" Brk (5) 


ersieht man, da die Multiplikation mit dem Tensor 


te = 8" 8x (6) 


den kovarianten MaBtensor des unverformten Koordinatensystems in den des verformten Ko- 
ordinatensystems tberfiihrt. Dieser Tensor! ist also kennzeichnend fiir den Formanderungs- | 
zustand. 


Aus verschiedenen Griinden ist es aber zweckmaBiger, den gemischtvariantenlogarithmischen | 
Formanderungstensor 


i 1 i 1 o; 
ek = 3 In (qi) = > In (8"" 8,4) (7) 


1 H. Tresca, Comptes rendus Acad. Sci. Paris 59 (1864) S. 754; Ann. du conservatoire des arts et métiers 
1865; Comptes rendus Acad. Sci. Paris 64 (1867) S. 809; Mém. prés. div. Savants 18 (1868) S. 733; Mém. prés. | 
div. Savants 20 (1872) S. 75 u. 281. 3 

2 M. Lévy, Comptes rendus Acad. Sci. Paris 70 (1870) S. 1323. 

R.v. Mises, Nachr. Wiss. Ges Géttingen 1913, S. 582; Z. angew. Math. Mech. 8 (1928) S. 161 

L, Prandtl, Nachr. d. Wiss. Ges. Géttingen 1920, S. 74; Z. angew. Math. Mech. 3 (1923) S. 401. 
Hi. Hencky, Z. angew. Math. Mech. 3 (1923) S. 241. 

R. Hill, The Mathematical Theory of Plasticity, Oxford 1950. 

W. Prager und P.G. Hodge, Theorie idealplastischer Kérper, Wien 1954. 

V.V. Sokolovskij, Theorie der Plastizitat, Berlin 1955. 

® DIN 1303, Schreibweise von Tensoren, Entw. 8. 57. 

10 Diese Schreibweise stimmt im wesentlichen mit der von A. Duschek und A. Hochrainer iiberein. 


1 A, Duschek und A. Hochrainer, Grundziige der Tensorrechnung in analytischer Darstellung, 3 Teile, | 
Wien 1950/55. 


12 4, E.Green und W. Zerna, Theoretical Elasticity, Oxford 1954. 
18 M. Lagally, Vorlesungen iiber Vektorrechnung, 5. Aufl., Leipzig: 1956. 
M4 Sein Zusammenhang mit dem Abbildungstensor und dem Dehnungstensor ist im Anhang 1 dargestellt. 
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einzufithren!. Schon H. Richter? wies auf die mit dieser Definition verbundenen Vorteile hin: 
1) Der logarithmische Formanderungstensor ist sinnvoll aufzuspalten in 


& =v + Bis (8) 
wobei der Formanderungsdeviator (Verzerrungstensor) 
i 1 1 esa 
Vk = EE — 3 ie E, (9) 
die Gestaltanderungen und der Kugeltensor 
i 1 i or 
fi => het (10) 


die Volumenanderungen verkérpert. 

2) Es gilt auch fiir endliche Formanderungen ein lineares additives Superpositionsgesetz. Des- 
halb kénnen z. B. die Gestaltanderungen in ihren elastischen und plastischen Anteil aufgespalten 
werden: 

: el pl. 
med ead (11) 

3) Der gemischt-variant definierte Formanderungstensor ist affin-invariant. 

H. Richter definiert den logarithmischen Formanderungstensor (7) ein wenig umstandlicher im 
ortsfesten Koordinatensystem. Die hier angegebene Definition geht in ihren Grundgedanken wohl 
auf H. Hencky* zuriick. Bei plastischen Formanderungen ist aber sonst anscheinend noch kein 
Gebrauch von dieser Definition gemacht worden, obwohl der Begriff ,,logarithmische Form- 
anderung“ in der Umformtechnik seit langem heimisch ist®. Er entspricht dem auf Hauptachsen 
transformierten logarithmischen Formanderungstensor. Entwickelt man ibrigens den logarith- 
mischen Formanderungstensor (7) in eine Reihe, so entspricht das erste Glied der in der modernen 
Elastizitatstheorie sonst iiblichen Definition des Formanderungstensors*. Der Zusammenhang des 


- log. Formanderungstensors mit dem Verschiebungsvektor ist im Anhang 2 dargestellt. Es ergeben 


aa 


cai alae aa 


sich nichtlineare Beziehungen. 

b) Invarianten des logm. Formanderungstensors. Fir die hier durchzufiihrenden Be- 
trachtungen ist es zweckmaBig, die Aufspaltung des logm. Formanderungstensors in Kugeltensor 
(Volumenanderungsanteil) und Deviator (Verzerrungstensor) beizubehalten. Fiir die Volumen- 
anderung ist nur die erste Invariante des Kugeltensors 


Bia pea ers 
en a (12) 
2 g V 


von Bedeutung. Hierbei ist g bzw. g die Determinante des kovarianten Maftensors (Gramsche 
Determinante) im verformten bzw. unverformten Zustand, die ein MaB fiir das Volumen V eines 
Kérperelementes ist. 

Beim Verzerrungstensor verschwindet die erste Invariante identisch 


@; =y'=0. (13) 


Die beiden anderen unabhangigen Invarianten werden hier in Anlehnung an Richter wie folgt 
definiert: 


On = yey! (14) 
und ; 
| Our = Ye Vi Yi (15) 
@,; kennzeichnet die GréBe der Gestaltanderung. Der Quotient 
i i 
praia Pa 16 
ea (16) 


1In (4;,) ist als Matrizenfunktion aufzufassen (vgl.Zurmiihl?) und nicht als der Logarithmus der Mafzahlen 


des Tensors. & 

2 R. Zurmiuhl, Matrizen, 2. Aufl., Berlin/Géttingen/Heidelberg 1958. 

3 H. Richter, Z. angew. Math. Mech. 29 (1949) S. 65. : 
_ 4 H. Hencky, Z. angew. Math. Mech. 5 (1925) S. 144, ferner C. Truesdell, Journ. Rat. Mech. and Analysis 1 
(1952) S. 125, 2 (1953) S. 593. : 

5 E. Siebel, Die Formgebung im bildsamen Zustand, Dusseldorf 1932. 

* Siehe FuBnote 12 von Seite 2. 

1* 
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ist also ein Ma fiir das Verhaltnis von Volumenanderung zur Gestaltanderung. Der Quotient 


Pin _ (7) | 


kennzeichnet dagegen die Art der Gestaltanderung. Die charakteristische Gleichung des Ver- |jf 
zerrungstensors 


i$ — 4 Gnd —> Ou = 0 (18) 


hat namlich nur reelle Lésungen fiir 0 < & < 1/6. Dabei gehért € = 0 zu einer reinen Gleitung, | 
€ = 1/6 zu einer einachsigen Verzerrung*. 


c) Geometrische Randbedingungen und Vertraglichkeitsbedingungen. Bei den | 
geometrischen Randbedingungen werden im allgemeinen nicht die Formanderungen, sondern die | 
Verschiebungen bzw. deren Richtungen gegeben sein: 

i = a(x tl nix, 4) oO (19) 


Die MaBzahlen des Verschiebungsvektors kann man auf das unverformte kérperfeste Koordinaten- | 
system beziehen (v') oder auch auf das verformte (v‘). Was man vorzuziehen hat, hangt vom | 
Einzelfall ab. 

Beim Ubergang von den Verschiebungen auf die Formanderungen, hat man zu beachten, daB 
sich in der Regel die Randbedingungen fiir die Formanderungen nicht ohne Willkir in einen 
elastischen und einen plastischen Anteil aufspalten lassen, zumal die Beziehungen zwischen Ver- 
schiebungen und Formanderungen nichtlinear sind. Probleme mit geometrischen Randbedingungen 
verlangen darum exakterweise nach einer gemeinsamen Behandlung der elastischen und plastischen 
Formanderungen. 

Ahnliche Uberlegungen gelten fiir die Vertraglichkeitsbedingungen, die sich aus der Forderung 
herleiten lassen, dafB der Riemannsche Krimmungstensor fiir das kérperfeste Koordinatensystem 
sowohl im unverformten wie im verformten Zustand verschwinden muB, da der Raum euklidisch 
ist und bleibt**. Die Vertraglichkeitsbedingungen betreffen den gesamten Formanderungstensor. 
Der elastische und der plastische Anteil werden fir sich allein in der Regel die Vertraglichkeits- 
bedingungen nicht erfiillen. Hierin liegt die Ursache der Eigenspannungen nach bleibenden Form- 
anderungen. Da die Beziehungen zwischen Kriimmungstensor und Formanderungstensor wiederum 
nichtlinear sind, ergeben sich auch hierbei schwer iibersehbare Verhaltnisse. 

d) Einige kritische Bemerkungen zur Definition des Formanderungstensors. Die 
Stetigkeitsforderungen an das Koordinatensystem zur Sicherung hinreichender Differentiierbarkeit 
bedeuten physikalisch, daB alle im unverformten Kérper auf einer stetigen Linie liegenden Punkte 
auch nach der Formanderung eine stetige Linie bilden miissen. Dem widerspricht die Erfahrung, 
daB alle bleibenden Formanderungen auf Versetzungen in den Gleitebenen der Kristalle beruhen. 
Um diesen Widerspruch zu beseitigen, sind Kontinuumstheorien der Versetzungen und Eigen- 
spannungen entwickelt worden, die von E. Kréner! ausfihrlich behandelt werden. Diese Theorien 
liefern einen verhaltnismafig einfachen Zusammenhang zwischen der Inkompatibilitat der auf 
Versetzungen beruhenden plastischen Formanderungen und den elastischen Restspannungen nach 
duBerer Entlastung. Dafiir sagen sie kaum etwas iiber den Formanderungsvorgang aus. AuBerdem 
sind sie meist auf kleine Formanderungen zugeschnitten. Zur Behandlung von Problemen der 
Umformtechnik sind diese Theorien daher im allgemeinen nicht geeignet. Immerhin sollte man sich 
bewubt bleiben, daB bei der hier vorgeschlagenen Beschreibung bleibender Formanderungen, die 
mit der Inkompatibilitat der plastischen bzw. elastischen F ormanderungen, sowie mit den geometri- 
schen Randbedingungen zusammenhangenden Fragen noch nicht befriedigend geklart sind. 


4, Spannungstensor. a) Definition des Spannungstensors. Die Spannungen miissen bei 
groBen Formanderungen stets auf den verformten Kérper bezogen werden. Ordnet man jeder Be- 
grenzungsflache eines aus dem verformten Kérper herqusgeschnittenen Elementes in bekannter 
Weise einen Vektor dF; so zu, da® sein Betrag der F lachengré8e entspricht und seine Richtung 
mit der nach auBen weisenden Flachennormalen n, zusammenfallt, so kann man fir die auf das 
Flachenelement wirkende Kraft schreiben 


IP Sordi (20) 


* Siehe FuBnote 3 von S. 3. 
** Siehe FuGBnote 12 von Seite 2. 


rock E, Kréner, Kontinuumstheorie der Versetzungen und Eigenspannungen, Berlin/Goéttingen/Heidelberg 
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D a7 ta L . . 
a ist oj, der auf das verformte Koordinatensystem bezogene Spannungstensor. Er ist sym- 
trisch, wenn am Kérperelement keine Momente angreifen. 
Der Spannungstensor ist aufzuspalten in 


att, 
= k= Tb Mh (21) 
. i pein i 
Up Cie at xe Oj, OF (22) 
der Spannungsdeviator und 
i 1 Deo, 
te = =z Of OF (23) 


der dem hydrostatischen Druck entsprechende Kugeltensor ist. 


5) Invarianten des Spannungstensors. Die Invarianten werden in gleicher Weise wie 
beim Formanderungstensor gebildet. Man erhalt also fiir den Kugeltensor die Invariante 


i= 0, = Gi, (24) 
die gleich dem dreifachen negativen hydrostatischen Druck ist: 
1 
Die Invarianten des Spannungsdeviators sind 
Yi=0, (26) 
Pu = Ti on Fy (27) 
Pur — <i on: z . (28) 
Yu kennzeichnet die Hohe des Deviator-Spannungs-Zustandes. Das Verhiltnis 
= ve 
wu =¢ (29) 


driickt aus, wie hoch die hydrostatische Beanspruchung gegeniiber dem Deviatorspannungs- 
Zustand ist. Der Quotient 
Pin 
Pit 


a, (o<o< 3) (0) 


gibt dagegen die Art des Spannungszustandes an, wobei ? = 0 zu einer reinen Schubbeanspruchung 
und # = 1/6 zu einem einachsigen Spannungszustand gehért. 

c) Gleichgewichts- und Randbedingungen. Fir die am Massenelement angreifenden 
Krafte gilt, wenn k; der Vektor der auf die Masseneinheit bezogenen Krafte ist und b; den Beschleu- 
nigungsvektor darstellt, die bekannte Gleichgewichtsbedingung? 


of, +0 (k;—b) =0. (31) 
Hier ist 9 die Dichte des verformten Werkstoffes. Sie laBt sich durch die Dichte im unverformten 


Kérper und durch die Determinante des kovarianten MaBtensors bzw. durch die erste Invariante 
des logm. Formanderungstensors ausdriicken: 


See 
re Shas (32) 
—— 0 & = 0 eB 
si =3 
Deshalb kann man fiir (31) auch schreiben 
oi en a a 
Eee ce! (33) 
bzw. a = 
Gl — o, Bi + hy —b =0, 
ist in Anhang 3 


ue oF ;, bezeichnet die kovariante Ableitung nach x*; vgl. Anhang 2. Die Ableitung von (31) 
L 


dargestellt. 
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wobei 
6, = t= Sh ee (34) 
g & 


die auf die Dichte bezogene Spannung ist. 
Fiihrt man in die Gleichgewichtsbedingung den hydrostatischen Druck und den Spannungs- |} 


deviator ein und bezieht auch diese auf die Dichte, indem man 


x pie P Xt rise Ti, 35 | 
-2,. u=4 (35) 
P Q : Q 
setzt, so erhalt man caat | 
tii — 4 Bl, — ple +p BL +h—h=0. (36) | 


Die Einfihrung der bezogenen Spannungen g' usw. ist zweckmaBig, weil hier alle GréBen auf | 


die Masseneinheit bezogen werden. Das liegt nahe, weil die Masse bei allen Formanderungen | 
konstant bleibt, das Volumen dagegen nicht. Die auf die Dichte bezogenen Spannungen haben | 
iibrigens die Dimension einer Energie je Masseneinheit. Das erhellt manche interessanten physikali- | 


schen Zusammenhange. 
Es ist angebracht, auch die Invarianten des Spannungstensors mit den auf die Dichte be- 


zogenen Spannungen zu bilden, also z. B. 


ie -s Py (37) 


einzufiihren. In gleicher Weise verfahrt man zweckmaBig auch mit anderen aus dem Spannungs- | 
tensor abzuleitenden Grifen. 

Die Randbedingungen lassen sich in der Form 

N; Oi = Sas) en, 6% my, = 0 (38) | 

angeben, wobei n; der Normal-Vektor der Oberflache und ¢‘*! der sogenannte ¢-Tensor ist. 


5. Allgemeine Betrachtungen von Zustandsanderungen. a) Energiegleichungen. Nach dem | 
allgemeinen Energiesatz gilt fiir abgeschlossene Systeme 
dW =6L+ 6Q. (39) 
Hierin ist + 
OW die Energieanderung des Systems, 
OL die mechanisch zugefiihrte Energie, 
6Q die als Warme zugefiihrte Energie. 


Die Energieanderung des Systems 1aBt sich aufgliedern in eine Anderung der kinetischen Energie 
OE und eine Anderung der inneren Energie 6U: 


“OW =6E + 0U. (40) 
Andererseits setzt sich nach den Grundgesetzen der Mechanik die zugefiihrte mechanische Arbeit 
oL teils in kinetische Energie 6E, teils in Formanderungsarbeit 54 um: 
6L =dE+ 6A. (41) 
Fiihrt man (40) und (41) in (39) ein, so folgt 
Soe | ‘ 


also 


6U = 6A + 60 


oder, wenn man bei gleichmaBiger Energiedichte die Energieanderungen auf die Masseneinheit 
bezieht, 


du = da + 0q. (43) 
Der Zuwachs an Formanderungsarbeit, bezogen auf die Masseneinheit, ist? 
da = 6; det , 
ae o 


* Die differentiellen Anderungen sind hier mit 6 gekennzeichnet, um anzudeuten, daB es sich nicht not-__ | 
ndig um Differentiale von ZustandsgréBen handelt. y 
® Ableitung in Anhang 4, 


(42) | 
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Bei ungleichmafiger Energieverteilung im Korper muB man von Gleichung (43) wieder auf (42) 
zuriickgehen. Man hat dann also beispielsweise zu bilden 


6A =f daodV. (45) 
V 


Die innere Energie u ist eine Zustandsfunktion des abgeschlossenen Systems. Ihre Anderung kann 
daher mit du bezeichnet werden (vgl. hierzu R. Becker*). Aus der allgemein giiltigen Energie- 
gleichung (43) ist aber nur eine Zustandsgleichung abzuleiten, wenn es gelingt, die rechte Seite 
dieser Gleichung durch ZustandsgréBen auszudriicken, von denen u eindeutig abhangt, so daB 
Zustandsanderungen unabhangig vom Wege sind. Das soll im folgenden kurz untersucht werden, 
wobei die Betrachtungen der Einfachheit halber auf die Masseneinheit bezogen werden. 


b) Reversible Zustandsanderungen. Bei reversiblen Zustandsanderungen ist 
6q = Tds. (46) 


T ist die absolute Temperatur und s die auf die Masseneinheit bezogene Entropie. Beide sind als 
ZustandsgréBen bekannt. Damit wird aus (43) 


du = 6; 6e; + Tds. (47) 
Beschrankt man sich nun auf adiabatische Zustandsinderungen (ds = 0) oder auf isotherme 
(T = konst), so kann man allgemeiner fiir (47) schreiben! 
dp — th oct, (48) 
wobei 
d du (adiabatische Zustandsanderung), 


ie d(u— T's) (isotherme Zustandsanderung) 


ein totales Differential einer Zustandsfunktion ist. Setzt man nun voraus, da y eindeutig vom 
logm. Formanderungstensor ¢j, abhange, ¢, also eine ZustandsgréBe sei, so folgt aus 


HEE aoe 
dy = oa dey, (49) 
die Integrabilitatsbedingung 
, nee (50) 
’ de}, 


oder bei Trennung von Volumendnderungs- und Gestaltanderungsarbeit 
x 0 xk 0 
—p=55 t= a (51) 
Die durch die Dichte 9 dividierten Spannungen leiten sich also im Raum der Formanderungen 
bei adiabatischen und isothermen Zustandsanderungen! — sofern die Formanderungen als Zu- 
standsgréBen anzusehen sind — von einem Potential ¢ (elastisches Potential) ab und sind deshalb 
eindeutige Funktionen der Formanderungen. Im einfachsten Falle ist bei isotropem Werkstoff 
das elastische Potential 


(Cal al 5 Scam ee 52 
Pe 2 Pu). (52) 
Hierin ist G der Schubmodul und m die Poissonsche Querzahl. B und ®,, sind die Invarianten 
des logm. Formanderungstensors. Der Zusammenhang zwischen logm. Formanderungen und 
Spannungen ist nach (52) und (51) 

: 2G m+1 


Z as 1D anicn ge eaen fae 
also é (53) 
—26m+1V, V 
a sae a area 
und 
‘ Xt 2G i 
ie ee 8 hh? 
} = 4, 
also (54) 


° 


i Vy, i 
Th =265y,- 


* Siehe FuGnote 2 von Seite 1. : : 
: Dascellie gilt noch allgemeiner fiir alle Zustandsanderungen, bei denen die Warmezufuhr eindeutig von. 


den Formanderungen abhangt. 
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Zu demselben Ergebnis kommt man, wenn man die durch die Dichte dividierten Spannungen |, 
als ZustandsgréBen einfiihrt. Dann ist 


Pee ty (35), 
dG; 1E 

pss es (56) | 

op 

5 é 
ee (57) | 

OFT | 

Aus 
ade 0 ii Xo Xi Xk ) 
°=46(8na1P +82), (58) | 
we m—2 x, A | 
=46 (3 =P + Vrs) 


folgen wieder die Beziehungen (53) bzw. (54). 
Bei kleinbleibenden Formanderungen gehen diese Beziehungen in die Hookeschen Gesetze tiber. | 
Die Erweiterung auf anisotrope Werkstoffe und allgemeinere Beziehungen zwischen Spannungen _ 
und Formanderungen bereitet keine grundsatzlichen Schwierigkeiten*. | 
Aus der Tatsache, daB das elastische Potential y eine skalare Funktion ist und sich diese Funk- 
tion sowohl durch den Spannungstensor wie auch durch den Formanderungstensor ausdriicken | 
laBt, folgt die Bedingung | 


oF ap? gyi = OP ap OP di. (59) 
/p Wa 


op op 
x OB ta avi, i 
d es dB und aia ap dy, (60) 
ap eth 


ist. Dies ergibt sich, wenn man einmal die Spannungsdnderungen als unabhangige Variable be- | 
trachtet und dabei die Formanderungsinkremente als abhangige Variable ansieht, das andere Mal || 
aber umgekehrt verfahrt. 

Aus diesen Beziehungen lassen sich die Extremal-Prinzipien der Elastizitatstheorie in all- 
gemeinster Form ableiten. So ist z. B. fiir isotrope Werkstoffe daraus zu folgern, da stets — also 
auch bei nichtlinearen Formanderungsgesetzen — 


Xs L ie 
dt, ~ dy, (61) | 
ist. Dies kann man auch so ausdriicken: Zustandsanderungen verlaufen so, daB die aktive Form- 
anderungsarbeit zu einem Maximum wird. 


c) Irreversible Zustandsanderungen. Bei allgemeinen Zustandsdnderungen, die auch 
irreversible Vorgange enthalten, kann 6g nicht mehr durch T'- ds ersetzt werden, da hierbei 


T ds > 6q ; (62) 
ist. Beschrankt man sich aber wiederum auf adiabatische Zustandsanderungen (dq: =-0) oder 


solche, bei denen die Warmezufuhr eindeutig von der Formanderungsarbeit anhangt (z. B. iso- 
therme), so behalt Gleichung (48) (mit df statt dy), also 


af = 8 bet | (48) 
ihre Giiltigkeit. Die Bedeutung des totalen Differentiales df andert sich allerdings im allgemeinen; — 
nur bei adiabatischen Zustandsanderungen ist wie vorher df = du. Der logm. Formanderungstensor 
bzw. der Spannungstensor kann aber jetzt nicht mehr fiir sich allein als ZustandsgréBe eingefihrt 
werden, da dann sofort wieder ein eindeutiger Zusammenhang zwischen Spannungen und Form- 
anderungen d. h. vollstandige Umkehrbarkeit der Vorgange zu folgern ware. 


Was man statt dessen als ZustandsgréfSen einzufiihren hat, und wie diese ZustandsgréBen von | 
Spannungen und Formanderungen abhangen, soll hier nicht weiter verfolgt werden. Es geniigt | 


* Siehe FuBnote 12 von Seite 2. 
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hier die Feststellung, da8 man bei allgemeinen Zustandsanderungen, die irreversible Vorginge 


enthalten, Spannungen und Formanderungen als unabhangige Variable anzusehen hat. Das 
bedeutet nicht, da8 Spannungstensor und logm. Formanderungstensor sich vollig unabhangig von- 
einander andern kénnen. Der Zusammenhang zwischen diesen beiden Tensoren wird aber nur 
eindeutig, wenn man gleichzeitig den Weg angibt, auf dem man die Spannungen oder die Form- 
anderungen dndert, wobei in manchen Fallen auch der zeitliche Verlauf anzugeben ist. 

Gleichung (48) behalt also nach diesen Ausfiihrungen fiir irreversible Zustandsanderungen nur 
den Charakter einer Energiegleichung. Diese kann man noch ein wenig umgestalten, indem man 
die Energieanderungen in ihren umkehrbaren Anteil dy (Anderung der potentiellen Energie) und 
in ihren nicht umkehrbaren Anteil dy (Dissipation) aufspaltet: 


df = dy + dy, ! 


x xi 63 
—— pdB + %, by;. (09) 


Wie sich die Gestaltanderungsarbeit auf potentielle Energie und auf Dissipation aufteilt, hangt 
von der Art der Formanderung ab. Die Volumenanderungsarbeit kann man dagegen in der Regel 
als vollstandig umkehrbar ansehen. Die potentielle Energie hangt eindeutig von den elastischen 
Formanderungen bzw. von den Spannungen ab. Sie ist also eine Zustandsfunktion der mechanischen 
ZustandsgréBen (Spannungen bzw. Formanderungen). 

Die dissipierte Energie ist dagegen keine eindeutige Funktion der Spannungs- bzw. Form- 
anderungen, sondern nur ihrer Differentiale. Sie ist im wesentlichen als Warmeenergie wieder- 
zufinden, teilweise tritt sie aber auch als latente Energie (z. B. in der Kaltverfestigung) auf. Fir 
die Dissipation gilt nach dem zweiten Hauptsatz der Thermodynamik 


oy =O. (64) 
Wo der zeitliche Verlauf eine Rolle spielt, geht man von 
) r . . ues. Xe 
ot 


aus, mit der Bedingung 
peo. (66) 


6. Forminderungsgesetze bei inelastischen Formanderungen. a) Ziahe Formanderungen. 
Zahe Formanderungen sind vom zeitlichen Ablauf, d.h. von der Formanderungsgeschwindigkeit 
abhangig. Die auf Newton zuriickgehende Beziehung 


te = 29k (67) 


-findet in der statistischen Thermodynamik durch die Onsager-Relationen* ihre Begriindung 


Verallgemeinerungen ergeben sich aus dem Ansatz** 
i oN 
ee (68) 
Oy; 
wobei N eine Funktion der logm. Gestaltanderungsgeschwindigkeit ist. 
Mit dem Ansatz (67) ergibt sich fiir die Dissipationsleistung eines Newton-Kérpers, bezogen 
auf die Masseneinheit, 


ee Qosi tk ; 
yp=2 Bo Vk Vi» : (69) 
= 2599. 
Allgemeiner ist 
| oan ta 2.20 
Vk 


Sind mit den zahen Formanderungen elastische verbunden, so ergeben sich verschiedene Még- 
lichkeiten des Zusammenwirkens. Die einfachsten Modellverstellungen findet man beim Se 
Kérper und beim Maxwell-Kérper, deren GesetzmaBigkeiten in der Ubersicht 1 angegeben sind**™*. 
Natiirlich lassen sich auch verwickeltere Arten des Zusammenwirkens darstellen. 


* Siehe FuBnoten 2 und 5 von Seite l. 
** Siehe FuBnote 5 von Seite 1. 
*** Siehe FuBnote 4 von Seite 1. 
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b) Idealplastische Formanderungen. Man spricht von ideal-plastischen Formanderungen, 
wenn ein Kérper sich plastisch deformiert d. h. in einen FlieBzustand gerat, sobald die Spannungen 
einen bestimmten Zustand (FlieBbedingung) erreichen, und wenn die FlieSbedingung nur die 
Spannungen betrifft und zeitlich unveranderlich ist. Der FlieBzustand dauert an, solange die 
Spannungen die FlieBbedingung erfiillen. 

Ist die FlieBbedingung 


so muf 
dee dn ao (72) 
OT}, 


sein, wenn der Korper im FlieBzustand bleiben soll. Nimmt man hierzu die — im Gegensatz zur 
Elastizitatstheorie stehende — Extremalbedingung, daB bei infinitesimalen Gestaltanderungen der 
Zuwachs der Formanderungsarbeit (aktiver plus passiver Anteil) zu einem Maximum werden 
solle, d.h. daB bei gegebenen dy* die Variation 


or}, dy; = 0 (73) 
sein soll, so folgt aus (72) und (73) das Formanderungsgesetz 
dy, = da (74) 
OT; 
baw. < 
+ nO, 


Hierbei ist dA baw. A ein beliebiger, im allgemeinen értlich und zeitlich veranderlicher, positiver 


Proportionalitatsfaktor. 
Der Ansatz (71) ist nicht auf isotropes Werkstoffverhalten beschrankt. Er ist ferner erweite- 
rungsfahig, indem man mehrere FlieBbedingungen ansetzt?: 


F(t) =0 (r=1...n). 7 (76) 


Von diesen FlieSbedingungen braucht immer nur eine erfiillt zu sein; es kénnen aber auch mehrere 
sein. Sind dies die Bedingungen ,F, so gilt 


Me er 
Paes ar? (Sm 1). ee) 


Die hier entwickelten Beziehungen werden sonst meist mit Hilfe des Begriffes vom plastischen 
Potential abgeleitet. Der so definierte ideal-plastische Kérper laBt sich als Grenzfall eines nicht- 
linearen Newton-Kérpers deuten*. Damit verlieren die sonst zur Begriindung der Theorie vom 
plastischen Potential erforderlichen Voraussetzungen ein wenig von ihrer Willkiir. Am haufigsten 
werden die v. Misessche FlieBbedingung, die dem Ansatz (71) entspricht, und die Trescasche FlieB- 
bedingung, die zu den verallgemeinerten FlieBbedingungen (76) gehért, verwendet. Hinsichtlich 
der letzteren, die in allgemeiner Form umstandlich anzugeben ist, sei auf Prager verwiesen. Die 
_ v. Misessche FlieBbedingung ist 


: Fai) =e — eo; (78) 
Zu ihr gehort das Formanderungsgesetz 
He =A te (79) 
Die Dissipationsleistung ist 
paar, 
ss (80) 
== Ach. 


Die Uberlagerung mit elastischen Formanderungen fithrt zu den Prandtl-Reufschen Gleichungen 
s. Tabelle 1). 


a oe a r veer » Probleme der Plastizitatstheorie, Basel-Stuttgart 1955; Proc. Inst. Mech. Eng. London 169 


* Siehe FuBnote 5 von Seite l. 
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Kin Hinweis sei hier noch angebracht: Die zeitlich und értlich unveranderliche FlieBbedingung 
(71) verlangt 


Fas (81) 


Diese Gleichung ergibt eine Bedingung, der die értlichen Ableitungen des Spannungsdeviators 
zusitzlich zu den Gleichgewichtsbedingungen (37) geniigen miissen. Der Spannungstensor ware 
damit iiberbestimmt, wenn sich seine MaBzahlen mit Hilfe der Vertraglichkeitsbedingungen ein- 
deutig auf drei reduzieren lieBen. Dies ist aber wegen des unbestimmten Faktors im Formanderungs- 
gesetz (75) baw. (77) nicht der Fall, und deshalb scheinen FlieBbedingungen und Gleichgewichts- 
bedingungen miteinander vertraglich zu sein. Diese Frage ist bisher aber noch nicht genau unter- 
sucht worden, obwohl eine Klarung winschenswert erscheinen mu; denn bekanntlich fiihren 
schon bei einfachen Problemen die obigen Ansatze manchmal zu Unstetigkeiten in der Spannungs- 
verteilung. 


c) Plastische Formanderungen mit Kaltverfestigung. Beriicksichtigt man bei plasti- 
schen Formanderungen die Kaltverfestigung, so wird die FlieBbedingung abhangig von der voraus- 
gegangenen Formanderung. Man hat deshalb die FlieBbedingung als Funktion des Spannungs- 
und des logm. Formanderungstensors anzusehen (vgl. Abschnitt 5 c)): 


Fein) 0 fase (82) 


Die Erhaltung des FlieBzustandes erfordert 


Fe 4 Se 04 (83) 
Ovi 


OTK, 


Stiitzt man_sich auf die gleichen Extremalprinzipien, wie sie fiir den elastischen Bereich gelten 


(vgl. Abschnitt 5 b)), so zerfallt Gleichung (83) in 
oF 
dz, = — 2e dyi. (84) 


= Xt 


Damit hat man das zur FlieBbedingung (82) gehérende Formanderungsgesetz gefunden. Bemer- 


_ kenswert ist, daB Gleichung (83) in Verbindung mit Gleichung (84) den Spannungsanderungen 


keine zusatzliche Bedingung auferlegt. Das mu8 auch gefordert werden; denn bei plastischen 
Formanderungen mit Kaltverfestigung sind die Beziehungen zwischen Spannungsanderungen und 
Formanderungen eindeutig im Gegensatz zu den Beziehungen bei ideal-plastischen Formanderungen. 

Fiir isotrope Kaltverfestigung kann man z. B., wie in der Theorie vom plastischen Potential 


es . 
ublich, von der Annahme ausgehen, daB die Invariante 1; des Spannungsdeviators eine Funktion 
der Gestaltanderungsarbeit sei 


Wee Ena) 0%. (85) 
Dieser Ansatz fiihrt zu dem Formanderungsgesetz 
xi 1 dk? i iomete xi 
dt, = oF da ok bzw. dy, = ah dt, . 186) 
da 


Das Inkrement des logm. Formanderungsdeviators ist danach dem Inkrement des Spannungs- 
deviators proportional. Eine Erweiterung des Ansatzes (85), die auch anisotrope Verfestigung 


- (Bauschinger-Effekt) beriicksichtigt, erhalt man z. B. mit 


Gap ny) r=, (87) 


wobei w und k2 im allgemeinen Falle noch Funktionen des Spannungs- und des Formanderungs- 


tensors sein kénnen. Dieser Ansatz, der bereits viele Anpassungsmdglichkeiten an das reale Werk- 
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Tabelle 1: Verh4 ae 


1 2 3 4 : 


6 
at Hooke Newton Kelvin Maxwell v. Mises” Prandtl-Rij 
a i i F F Jastisch || 
Formanderung elastisch eek (mE bey ties isco ideaiplastl ie 
einachsiges 4 p y p j p |? p p : 
mechanisches jr H = EH H EW AS ee 
Modell | 
FlieBbedingung 
Volumen-| * 2G m+1 B he wie lL wie 1 
anderg. 36 m—2 
Form- vice 10.6 ae 0%. 
anderungs- | ti = —— yf | dh =a th | 
gesetze Gestalt- ie ae A, Ka oe g 26 : 
2 t= —— yh t= —— vk 9 - 
anderg. fa) Q er Pa @ x; 
Be Bias | 
A a na le | | 
30 m—2 | . ie 
Volumen- 5 eae eV te wie 1 SSeS wie I ji} 
anderg. b) 30 m—2 x, | } 
Potentielle 4G m+1 
Energie : i 
f G i wk 1 i 
a) = yp vk wie la | 
Gestalt- Q | 
anderg. : j 
b) fe aj te wie 1b wie 1 bil) 
OSE ecole : eas We 
a) —— vi v4 wie 2a Ati te =Ak 
Dissipationsleistung e re. wie 
7 2 = xt 
%, b) 5 % wie 2b 
n 


Die Unterlage des Kérpers ist in der Weise nachgiebig zu denken, da er mit zunehmendem Weg (ohne Riicksicht auf die Richtung) tiefer einsinkt 


) 


stoffverhalten bietet, ist vom Verfasser bereits an anderer Stelle! behandelt. Das zugehérige Form- 
anderungsgesetz ist 


ou ak? : 
2u(i—n yt) +2 55 vi (ie — pw yt) + avi 
i 3 
dt = ae 5 du agree a ok2 dy, (88) 
2(%—pyt)—2 xi Yn A or XE 
OT}, OT 


Der Kreis der Méglichkeiten ist damit jedoch keineswegs erschépft. Auch eine Uberlagerung — 


elastischer Kigenschaften in der iiblichen Weise bereitet keine grundsatzlichen Schwierigkeiten 
(s. Tabelle 1). Desgleichen lassen sich zahe Formanderungen tiberlagern. 


d) Kritische Stellungnahme zu den Gesetzen plastischer Formanderungen. Die 
hier vertretene Auffassung, daB bei plastischen Formanderungen mit Kaltverfestigung Spannungs- 
tensor und logm. Formanderungstensor (nicht ihre Differentiale) als unabhangige Variable anzu- 


1 Th, Lehmann, Z. angew. Math. Mech. 38 (1958) S. 295. 
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einfacher Modellkérper 


ee ee 
ls 8 9 ; 10 11 
plastisch mit is = = 
5 plastisch mit anisot: A x . 5 . 
cating | MERRIE” | reticent) | seek pi | nt 
1) 2) 1) 2) 1) 2) F 1) 2) 
ld Z a j ‘ Pp ; y 
oes #3 (ti —u yi) Z 
tit —k(a) = 0 dt F(z; yj) = 0 wie 9 deacon Res 
x (7 — yt) —k =—0 5 Veo VE] = 
= 
wie 1 wie 1 
gaa — t= (g6+4 
ea igs ne (ae, = (ac = i i 
Oy Beer eT Been kK oe k x 
da Vig ha said os x (4-F vi 
= A(a) ti x. Neos see Ati vi) x 2 Ui 
ce = + Ati; yj) ti hm a(t 8 VR vi 
wie 1 wie | 
a prone 27m | 
AG AE 6 i) 
wieaah x (4 i) 
t 6 t 
EA Aaj vs Att TH —¢ 
A = Xa) A = konst. 7 A(zis vi) wie 9 


* Die in der nachgeschalteten Feder gespeicherte Energie ist nicht mehr als potentielle Energie anzusprechen, sondern richtiger als latente 
Energie, die nur bei Belastungsumkehr wieder frei wird. 


sehen sind, steht im Gegensatz zu den bisher tblichen theoretischen Ansatzen*. Diese gehen von 
der FlieBbedingung fiir ideal-plastische Formanderung aus, also z. B. von der Gleichung (79); 


man betrachtet lediglich die invariante GréBe k2 als abhangig von der Formanderungsgeschichte. 
Der Begriff vom plastischen Potential wird beibehalten und 


dy = aie (89) 


OT; 


-gesetzt, wobei F und A Funktionen sind, die nur vom Spannungstensor abhangen, wahrend ie 
eine Funktion der Gestaltanderungsarbeit a ist. Die Bedingung fiir die Erhaltung des F'lieB- 
-gustandes fiihrt dann schlieBlich unter Elimination der Funktion h zu dem Formanderungsgesetz 


1 oF 


* Siehe FuBnoten 6 von Seite 2 und 1 von Seite 10. 
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x 
wobei F' als homogene Funktion n-ten Grades angenommen ist. Speziell wird mit F=Y¥y 


dye oe age (91) 
YI : 


Diesen Gedankengiingen liegt jedoch eine falsche Beurteilung des physikalischen Sachverhaltes 
zugrunde. Hill* geht beispielsweise davon aus, da zu der Bestimmung der 6 Formanderungs- 
inkremente 6 Gleichungen bendtigt werden. Dazu rechnet er 3 Gleichgewichtsbedingungen, 1 FlieB- 
bedingung, 1 Bedingung iiber die GréBe der Formanderungsinkremente als Funktion der Spannungs- 
inkremente, 1 Inkompressibilitatsbedingung. Hierzu treten die allgemeinen Formanderungs- 
gesetze (74) der Theorie vom plastischen Potential, die die Verbindung zwischen Spannungen und 
Formanderungen vermitteln. 

Demgegeniiber ist festzustellen, daB zur Bestimmung der Formanderungsinkremente nur 
3 Gleichungen bendtigt werden; denn von den 6 Formanderungskomponenten lassen sich stets 3 
mit Hilfe der Kompatibilitatsbedingungen eliminieren. Darin kommt die physikalische Tatsache 
zum Ausdruck, daB sich der Formanderungstensor auf die 3 Verschiebungskomponenten zurtick- 
fiihren lat. Zu ihrer Bestimmung reichen! die 3 Gleichgewichtsbedingungen (34) in Verbindung 
mit den Forminderungsgesetzen (84) aus. Greift man bei den Gleichgewichtsbedingungen auf die 
Gleichung (37) zuriick mit der Aufspaltung in Spannungsdeviator und Kugeltensor, so tritt hierzu 


noch die Bedingung 


weil die Gleichgewichtsbedingungen in diesem Falle eine Unbekannte (p) mehr enthalten. Die 
Inkompressibilitatsbedingung ist durch die Formanderungsgesetze (84), die sich aus der Flie- 
bedingung (82) ableiten, eliminiert. Die FlieBbedingung ist bei dieser Betrachtung, wie schon ein- 
mal hervorgehoben wurde, keine zusadtzliche Bedingung; sie bestimmt nur die Formanderungs- 
gesetze. Insofern ist sie allerdings von zentraler Bedeutung. Wieweit die FlieBbedingung dadurch 
eingeschrankt ist, daB die mit Hilfe der zugehérigen Formanderungsgesetze aus den Gleichgewichts- 
bedingungen ermittelten Formanderungen den Kompatibilitatsbedingungen geniigen miissen, ist 
bisher nicht untersucht. Die Klarung dieser Frage vermag sicher noch interessante Einblicke zu 
vermitteln. 

Fir die bisherige Betrachtungsweise lieBe sich noch anfithren, da experimentelle Versuche zur 
Klarung der Formanderungsgesetze im allgemeinen eine befriedigende Ubereinstimmung mit der 
Theorie vom plastischen Potential ergaben”*. Bei manchen Versuchen ergaben sich allerdings 
auch gewisse systematische Abweichungen, wenn auch von geringer Gréfe*. Theoretische Unter- 
suchungen fithrten dagegen auch schon zu Widerspriichen®. Fir die neue Betrachtungsweise liegen 
noch keine Ergebnisse theoretischer oder experimenteller Untersuchungen vor. Die Ansatze lassen 
aber so viel Freiheiten offen, daB es méglich sein sollte, die Theorie weitgehend dem realen Werk- 


_stoffverhalten anzupassen. Fir Belastungsfalle, bei denen dtj. stets proportional Th ist, fuhren 
iibrigens beide Betrachtungsweisen zu den gleichen Ergebnissen®. 


Auf einen Punkt sei abschliefSiend noch einmal hingewiesen. Die Gleichung (84), die zur Her- 
leitung des Formanderungsgesetzes benutzt wird, enthalt noch eine gewisse physikalische Willkiir, 
auch wenn sie sich an die Gesetze elastischer Formanderungen anschlieBt. Es ware denkbar, diese 
Gleichung fallen zu lassen. Die Forderung, da8 FlieBbedingung und Formanderungsgesetze so 


formuliert werden miissen, da® Gleichung (83) stets erfiillt ist und sich kein Widerspruch zu den_ |} 


Gleichgewichtsbedingungen ergeben kann, bleibt aber davon unberihrt. 


7. Beanspruchung der Werkstoffe. Bei der Kaltumformung handelt es sich vorwiegend um 


elastisch-plastische Formanderungen mit Kaltverfestigung, die im allgemeinen stark anisotrop ist. | 


Bei der Warmumformung machen sich dagegen bereits Zahigkeitseinfliisse starker bemerkbar. 


‘ 


*Siehe S. 33 der in FuBnote 6 von Seite 2 zitierten Arbeit. 

1 Wegen der Eindeutigkeit der Beziehungen zwischen Spannungsanderungen und Formanderungen bei 
plastischen Formanderungen mit Kaltverfestigung. 

® W. Sautter, A. Kochendérfer und U. Dehlinger, Z. Metallkunde 44 (1953) S. 442 und 553. 

° U. Dehlinger, J. Diehl und J. MeiBner, Z. Naturforschung 1la (1956) S. 37. 

* J. M. M. Morrison und W. M. Shepherd, Proc. Inst. Mech. Eng. London 163 (1950) S. 1. 

5 H. Lippmann, Z. angew. Math. Mech. 38 (1958) S. 297. 

* Das plastische Potential bei ideal-plastischem Verhalten la8t sich ferner auch als Grenzfall der neuen 
Betrachtungsweise deuten, steht also nicht im Widerspruch dazu. 
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Rein viskose Stoffe zeigen keine Erschépfung des Formanderungsvermégens. In der Warm- 
umformung driickt sich dies im Rekristallisationsvermégen der Werkstoffe aus. Die Rekristalli- 
sationsgeschwindigkeit hat allerdings eine endliche GréBe und ist zudem stark temperaturabhangig. 
Daher kann man die Warmumformung nicht als rein viskos ansehen. Zur Kennzeichnung der 
Beanspruchung geniigen deshalb nicht allein die Invarianten des Spannungstensors bzw. des 
Tensors der logm. Formanderungsgeschwindigkeit. Noch starker ist es bei der Kaltumformung aus- 
gepragt, daB die GréBe der Beanspruchung nicht allein durch einen Tensor gekennzeichnet werden 
kann, wie schon die thermodynamische Betrachtung der ZustandsgréBen ergab. Zwar kann man 
den Ubergang vom elastischen zum plastischen Bereich allein durch den Spannungstensor festlegen; 


x 
bei isotropen Werkstoffen ist es bekanntlich vor allem die Invariante Wit, von der dieser Ubergang 
abhangt. Die Erschépfung des Formanderungsvermégens hangt dagegen auBerdem! vom hydro- 


statischen Druck p und vor allem auch von der voraufgegangenen Formanderung, also vom logm. 
Formanderungstensor, von der Gestaltanderungsarbeit usw. ab. Eindeutige ZustandsgréBen zur 
Kennzeichnung einer Beanspruchung sind noch nicht ermittelt. So mu8 man sich vorerst darauf 
beschranken, fiir bestimmte Belastungsvorginge — fiir diese sind die Formanderungen eine ein- 
deutige Funktion der Spannungen und umgekehrt — die Formanderungsgrenzen anzugeben. Zu 
ihrer Kennzeichnung geniigt dann die Angabe des Formanderungszustandes bzw. des Spannungs- 
zustandes. Dabei bedient man sich vorteilhaft der Invarianten und der aus ihnen abgeleiteten 
Verhaltniszahlen (s. 232 und 242). Bei anisotropen Werkstoffen gehen allerdings in die kenn- 
zeichnenden Invarianten auch gewisse Kérpereigenschaften ein. 


Anhang 1 


Abbildungstensor, Dehnungstensor und logm. Férmanderungstensor. Die Umgebung eines K 6rper- 
punktes mit den kérperfesten Koordinaten x‘ bildet sich auf ein ortsfestes, cartesisches Koordi- 
natensystem y’ im unverformten Zustand (gekennzeichnet durch eine tibergesetzte 0) mit Hilfe 
der Transformationsformeln 


dy' = ¢, dx* (1.1) 
ab, deren Umkehrung 

dx = ¢; dy' (1.2) 
ist. Im verformten Zustand lauten die Transformationsformeln 

dy? = ci dx (1.3) 
bzw. 

dx* — cf dy’ . | (1.4) 


Die Transformationsmatrizen sind im allgemeinen Funktionen des Ortes und der Zeit, sie sind 
keine Tensoren, weil sie sich auf zwei verschiedene Koordinatensysteme beziehen. 
Bildet man dagegen im ortsfesten Koordinatensystem den verformten Kérper auf den unver- 
formten ab mit Hilfe von 
dy = cj, dx* , 
iP 2h 70 
= CC, dy’ ; (1.5) 
= Ai, dy*, 


so erkennt man, daB diese Abbildung durch den gemischt-varianten Abbildungstensor Ai, ge- 


Jeistet wird, der auf das ortsfeste Koordinatensystem bezogen und im allgemeinen unsymmetrisch 


ist. Er 1a8t sich multiplikativ aufspalten in einen symmetrischen Dehnungstensor Di" baw. Dj, 
oder D,,, der eine drehungsfreie Formanderung vermittelt, und in einen Versor R;, bzw. R', oder 
Ri*, der eine reine Drehung bedeutet. Bec 45 ; 
Die Aufspaltung ist in verschiedener Weise méglich, je nachdem, welche Operation man sich 
zuerst ausgefihrt denkt und auf welche Grundvektoren man die Tensoren bezieht. Formal ergeben 
sich dabei zwei verschiedene Zerlegungen, von denen jede wiederum zwei verschiedene Deutungen 
zulaBt: : 
Ag Ds Rk pp Rees (1.6) 
Ai, = Ri, DP = Ri” Dia, - (1.7) 


1 E. Siebel und H. BeiBwanger, Tiefziehen, Miinchen 1955. 
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Die Schreibweise (1.6) laBt folgende Auslegungen zu: ; 
(1.6 a) Man hat sich zuerst die Drehung ausgefiihrt zu denken und dann die Formanderung. 
Beide Tensoren sind dabei auf das ortsfeste Koordinatensystem bezogen. 


(1.6b) Man hat sich zuerst die Formanderung ausgefiihrt zu denken, wobei das ortsfeste 
Koordinatensystem mit transformiert wird (Basistransformation). Der Versor der anschlieBenden 
Drehung ist auf das transformierte Koordinatensystem bezogen. 


Fiir die Schreibweise (1.7) erhalt man folgende Deutungen: 


(1.7 a) Man hat sich zuerst die Formanderung ausgefiihrt zu denken und dann die Drehung. 
Beide Tensoren sind dabei auf das ortsfeste Koordinatensystem bezogen. 


(1.7b) Man hat sich zuerst die Drehung ausgefiihrt zu denken, wobei das ortsfeste Koordinaten- 
system mitgedreht wird (Basistransformation). Der Tensor der anschlieBenden Formanderung ist 
auf das gedrehte Koordinatensystem bezogen. 


Hier wird den beiden Aufspaltungsméglichkeiten jeweils die erste Bedeutung beigemessen. 
Man kommt aber zu demselben Ziel, wenn man jeweils von der zweiten Bedeutung ausgeht. 


Aus (1.6) eliminiert man den Versor, indem man A’, von rechts mit dem konjugierten Tensor 
Aj} multipliziert. Das ergibt nach den Regeln fiir die Bildung konjugierter Tensorprodukte (Ver- 
tauschung der Reihenfolge und Ubergang zu den konjugierten Tensoren) : 


Q! = Ai, Ay’ = Dan, Re De, 
(1.8) 


ray po k 
ce il (RAT Set Da 
Da der Versor im ortsfesten cartesischen Koordinatensystem ein orthogonaler Tensor ist, fiir den 


ne Rie an On Om baw. Tags (a ai Sah One (1.9) 


ist, folgt also 


Qi* = Di df, 67, DE = Di, Dis 
(1.10) 


ayy k — Di pk 
= D*’ 6m Oms D® ry DD 
Q'* ist also das Quadrat des symmetrischen Dehnungstensors, bezogen auf das ortsfeste Koordinaten- 


system, fiir den Fall, da8 die Drehung vor der Formanderung ausgefiihrt zu denken ist. Setzt 
man fiir A‘,, baw. A;* die Werte aus (1.5) ein, so findet man 


tk i ok itr = k 
ae Ae Aig ie) Cape 0. 
e : (1.11) 


SoG ere ck 
= ¢ g"* c,. 
Transformiert man nun diesen Tensor auf das verformte kérperfeste Koordinatensystem, was bei 


der hier zugrunde gelegten Bedeutung des Dehnungstensors nahe liegt, so ergibt sich nach den 
Transformationsregeln ; 


(1.12) 


pe 2 WUT 3 eee nk 
q =¢,Q Ch = Cm & ast 


= Oi gr ok = BF. 


. Der gemischt-variante quadratische Dehnungstensor, bezogen auf das verformte kérperfeste Koor- 
dinatensystem, ist also 


qe = 9" Br = 8 Brn (1.13) 
Den logarithmischen Formanderungstensor erhalt man aus ihm durch die Beziehung 


ek, = In Vg) => In (qi), (1.14) 


wobei der Logarithmus und die Wurzel als Tensorfunktionen aufzufassen sind. 
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Geht man von der Darstellung (1.7) aus, so erhalt man durch Multiplikation mit dem konjugierten 
Tensor von links einerseits 


0: a A; Avs a 
ir Rom Re Dies (1.15) 


also wieder das Quadrat des Dehnungstensors bezogen auf das ortsfeste System, freilich jetzt fir 
den Fall, da die Formanderung vor der Drehung ausgefiihrt zu denken ist. Andererseits ergibt 
sich durch Einsetzen 

~ Se pi eee SS 2 

Qn =Ai A= Cre, c. Ck 
(1.16) 


° 
=> =) 


—s 
= Ci 8rs Ck 


und nach Transformation auf das unverformte kérperfeste Koordinatensystem, was bei der hier 
zugrunde gelegten Bedeutung des Dehnungstensors allein sinnvoll ist, 


; : (1.17) 
= 0; 8,5 Ok = Bik 
Daraus folgt wieder 
Gi = 2" Gen = 8°" 8,4 - (1.18) 
Es ist also 
CRESU He (1.19) 


_d.h. es ergibt sich aus beiden Betrachtungsweisen formal das gleiche Ergebnis. Das gilt allerdings 
nur, wenn man die Tensoren gemischt-variant schreibt. Damit ist die Bevorzugung dieser Schreib- 
weise bereits hinreichend begriindet. 

Geht man von den beiden anderen Deutungen aus (1.6 b) und (1.7), die die Schreibweisen 
(1.6) und (1.7) noch zulassen, so wird man ebenfalls zu den formal gleichen Ergebnissen (1.13) 
bzw. (1.18) geftihrt. Der logarithmische Formanderungstensor (1.14) laBt also vier verschiedene 
physikalisch sinnvolle Deutungen zu, wobei er teils auf das unverformte, teils auf das verformte 
kérperfeste Koordinatensystem (teils vor der Drehung, teils nach der Drehung) bezogen zu denken 
ist. Hierin d4ufert sich seine affine Invarianz. 

Die Bedeutung des Tensors qi, ergibt sich ferner besonders daraus, da8 er durch Multiplikation 
den kovarianten MaBtensor des unverformten kérperfesten Koordinatensystems in den des ver- 
formten iiberfiihrt 


Bim Ue = Bim 8" Brk = Bik: (1.20) 
Der inverse Tensor 
(qi)? = 8" Brn = Bur 8" (1.21) 
uberfiihrt hingegen den kontravarianten Mafitensor 
eGie =f" 55:8" = 8" (1.22) 
Anhang 2 


Zusammenhang zwischen F ormanderungstensor und Verschiebungsvektor. Im ortsfesten carte- 
sischen Koordinatensystem hat der Verschiebungsvektor u‘ die MaBzahlen 


7 ui = yi — yi, (2.1) 
Es ist also 
ice SSL ae 
OX}, axk Oxk ? 
q — Ch — rot Py (2.2) 
oder ; a dui 
Ch = & Tae 


= nieur-Archi 
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Multipliziert man die Matrix cj von links mit der konjugierten Matrix, so erhalt man 


ou™ a, du™ du™ du™ 
L om : 
it Ggh ok Gale eat pak 


also 
7) 


m m 
DN ata 3 4 am du™ L am du™ ; du’ ou 
Sik — Sik i Oxk k Oxi | Oxi dxk 


(2.3) 


Transformiert man u‘ auf das unverformte kérperfeste Koordinatensystem, so folgt 
a = ci uk oder “=e, v". (2.4) 
Mit Hilfe des absoluten Differential-Kalkiils* findet man sodann 


*) am 
Ge 2 Gls dcm ~, 


oxt = Oxi gage 2 


(2.5) 


ye Bee oye 
salt sax ” axt ‘ 
* (2.6) 
ovr Sr ~s 


Ij, sind die auf das unverformte Koordinatensystem bezogenen Christoffelschen Symbole zweiter 
Art, die sich beispielsweise in folgender Art schreiben lassen: 


(2.7) 


Lorn (in 1 hen Bie) 
a Oxs Oxt Oxk 


Setzt man (2.5) in (2.3) ein, so erhalt man 


° ° ° bod ° ° ~ } 
Be Sie ee ee rs ix 
° ° ~ ° | ° ~o) ~T\ 
= Sik “Sir pot Bee 8s, et v" wile (2.8) 


ete eh NS Rak rahe Be 
= Bin + Ml, +), + Od, 2], 


und schlieBlich 
qi = 8°" Sn = Of +O], + 0,|' +0] Uh (2.9) - | 


Mit qj, kann man nun den logarithmischen Formanderungstensor 


5 1 . 
ei = + In gi) (2.10) 


bilden und hat damit den gesuchten Zusammenhang zwischen Formanderungstensor und Ver- 
schiebungsvektor. 


Man kann den Verschiebungsvektor auch auf das verformte kérperfeste Koordinatensystem | 
bezichen. In diesem Fall lést man Gleichung (2.2) zweckmaBig nach cj, auf und multipliziert dann _ | 
wiederum mit der konjugierten Matrix von links. Nach Transformation des Verschiebungsvektors — | 
auf das verformte Koordinatensystem erhalt man zunichst . 
| 
| 
| 


rk = Sin — Ui, — Yn; ae sels Dial : (2.11) 


* Siehe FuBnoten 12 und 13 von Seite 2. 
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und schlieBlich 


Sgn Ge) — 6, — 0 yl |! ol (2.12) 
Den logarithmischen Formanderungstensor erhalt man daraus durch die Beziehung 
i : i 1 
Sh 9. In (Gj) = — > In (qi). (2.13) 
Anhang 3 


Gleichgewichtsbedingung. Es wird das Kraftegleichgewicht an einem Parallelepiped des ver- 
formten Kérpers betrachtet, das durch die Ebenen x‘ = konst. und xi + dx‘ = konst. gebildet 
wird. Die durch die Ebenen x‘ = konst. gebildeten Begrenzungsflachen werden durch die Vektoren 


dF, = — /g dx™ dx”, (m == n =+ i) 
aeea Vy" (3.1) 
ame 
reprasentiert, wobei 
dV =)/g dx! dx? dx3 (3.2) 


das Volumen des Parallelepipeds und g die Determinante des kovarianten Maftensors im Punkte 
: I aos 2 3 
ot = dx‘ ist. Zu den Ebenen x‘ + dx‘ = konst. gehéren die Vektoren 


dF, = + Vg dx™ dx", 


Sen ei ides 
te (m + 1) (3.3) 


Fiir die an diesem Volumenelement angreifenden Krafte gilt nun 
o, dF; + (0, + 0; , dx’) dF; + odVk,—odVb,=0 (i summiert) . (3.4) 


Hierin ist @ die Dichte des Kérpers, k, der Vektor der auf die Masseneinheit bezogenen Krafte 
und 6, der Beschleunigungsvektor. 


Setzt man dF, und dF; ein, so folgt 
; dV - dV : ) 
=O, ae Cbg 1 Oey ot OV (hy —b,) = 0 | 


oi]; +0 (& — b,) = 0. 


(3.5) 


oder 


Das ist die bekannte Gleichgewichtsbedingung. 
Bildet man das Moment der Krafte um den Mittelpunkt des Parallelepipeds, so erhalt man 
unter Vernachlassigung der Glieder, die von héherer Ordnung klein sind, die Bedingung 


oe, 5°" oe . = dxk = 0 (k pect (3.6) 
oder A 
e410 = 0, 
3.7) 
2 nao : 


Daraus folgt also die Symmetrie des Spannungstensors. 


Anhang 4 


Formanderungsarbeit. Unterwirft man die Volumenelemente eines Kérpers, der sich im Gleich- 
‘gewicht befindet (b; = 0), virtuellen Verschiebungen dv‘ und integriert tiber das gesamte Kérper- 
volumen V, so erhalt man 
f(cl. év' + gk, dv‘) dV = 0 
ah 4.1) 
d . ( 
a f (orl dv, + @ ki dv;) dV =0. 
eae 
9* 
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Der Vektor der virtuellen Verschiebung ist hierbei wie oj, bzw. o'* auf das verformte Koordinaten- 
system zu beziehen. 
. . . . ik 
Daraus folgt zunachst, wenn man die kontravariante Schreibweise des Spannungstensors 0° 
benutzt, die hier etwas zweckmaBiger ist, 


Wed ed ae ie dv), dV + fe ki dv, dV = 0. (4.2) 
Vig V ie =a ¥: 


Formt man das erste Integral mit Hilfe des Gaufschen Integralsatzes in ein Oberflachenintegral 
um, so erhalt man 


fo du; dF’, + fo ki du,dV = fot ov,, dv. (4.3) 
F v Vv 


Der erste Term der linken Seite stellt die Arbeit der Oberflachenkrafte, der zweite die Arbeit der 
Volumenkrifte dar. Die rechte Seite ergibt die Formanderungsarbeit. 


Wegen der Symmetrie des Spannungstensors kann man fiir die rechte Seite auch schreiben 


[ot dei av =z fot (dui, + de.) dV. (4.4) 
74 


Fiir den Zusammenhang zwischen Verschiebungsvektor und MaBtensor gilt nach Anhang 2, Glei- 
chung (2.8) allgemein 


Sik — Sik = Asix = v;, si v;, ; iz ,| ; vl, ’ (4.5) 


wobei auf der rechten Seite alles auf das unverformte kérperfeste Koordinatensystem bezogen ist. 
Beachtet man, da hier nur virtuelle Verschiebungen betrachtet werden, bei denen das verformte 
Koordinatensystem als Bezugsbasis dient, so ergibt sich 


Ogin = Ov;|, + Ov,|, + dv, , dv", . (4.6) — 


Das letzte Glied der rechten Seite ist bei infinitesimalen Verschiebungen vernachlassigbar klein 
Ks wird also 


Ogi, = Ov], + Oy; , (4.7) 


und fiir die Formanderungsarbeit erhalt man 


Ib : 
OA = x | ods dV. (4.8) 
Vv 


In dieser Form wird die Formanderungsarbeit iiblicherweise in der modernen Elastizitatstheorie 
angegeben*. 


Geht man nun wieder zum gemischt-varianten Spannungstensor iiber, so erhalt man zunachst 


1 i 
- bd = > [ ot, gmk by, dV (4.9) 
y 


oder nach Erweiterung mit dem MaBtensor des unverformten Koordinatensystems 


1 4 ° ° 
04. i Om (8 Srn) O(8”"* 8,;) dV. (4.10) 


v 
Benutzt man die in Anhang 1 abgeleiteten Definitionen fiir den quadratischen Dehnungstensor 
gg =s"'g; und (q3) = g""g.,, (4.11) | 


* Siehe FuBnote 12 von Seite 2. 


aim = : d ct 
ty me 3 
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80 geht (4.10) aber Ate, 2 aga ee 


I| 


es ce ad x | om (any 697 AY, | z 


SS Sie : = | om > dlm(ar) dV, (4.12) 


ey 


ae gear 3 = | oi, de" dV, | 


bei ej der logarithmische Formanderungstensor ist. mer: 
Die Formanderungsarbeit eines Elementes, bezogen auf die Masseneinheit ist schlieBlich 


pe eee i a == oi, de” = oi, Oe? (4.13) 
aettue paamenete sl ai. S. 4.14) 
= : Om == 0 OF . = ( Ms ) 
: = __ (Eingegangen am 9. April 1959. eran) : 


Verfassers: Dozent 


Dr.-Ing. Theodor Lehmann, Hannover, PodbielskistraBe 93. 
ee ee eet. See 4 i iss * 
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Uber die Biegung einer unendlich erstreckten orthotropen Platte 
auf elastischer Unterlage 


Von S. Woinowsky-Krieger 


1. Grundgleichungen. Es liege eine Platte konstanter Dicke vor, die in bezug auf ihre elastischen 
Eigenschaften iiberall zwei zueinander und zu ihrer Mittelebene senkrechte Symmetrieebenen be- 
sitzen mége. Sind x = konst. und y = konst. diese Ebenen, so lautet die Differentialgleichung fiir 
die Einsenkung einer solchen orthotropen elastisch gebetteten Platte 
o*w 
Ox! 


o*w d*w 
Vw = D, aa + 2H aaa + Dy Ga = UY) — PY)» (1) 
Hierin ist q die Intensitat der gegebenen Auflast und p die Intensitat der Bodenreaktion, wahrend D,, 
H, DD; Plattensteifigkeiten bezeichnen, deren Bedeutung wir als bekannt voraussetzen wollen.* 

Die Eigenschaften der Unterlage lassen sich in einer ziemlich allgemeinen Weise durch die An- 
nahme erfassen, daB eine an der Stelle (£, 7) der Bodenoberflache aufgebrachte Punktlast P = 1 
an einer Stelle (x, y) dieser Oberflache eine Einsenkung K(s) erzeugt, wobei s = V(x — &)? + (y—n)? 
und K eine von den Eigenschaften der Unterlage abhangige Funktion von s bezeichnen soll. Der 

‘ gegenseitige Druck p(é, 7) zwischen Platte und Unterlage erzeugt somit an der Stelle (x, y) eine 
Gesamteinsenkung 


wins) I plea) EAs) dadg— late, ah Vales al heh ae ar, (2) 


wobei das Doppelintegral sich iiber die gesamte Begrenzungsebene, d. h. von £ = —oo, 7 = —oo 
bis € = + 00, 7 = + on, erstreckt. 

Unter der Annahme einer zum Koordinatenkreuz doppeltsymmetrischen Belastung q setzen 
wir zur Auflésung der Gleichung (2) 


w(x, y) eA B) cosa x cos B y dx dB (3) 


_ und, sofern das nachstehende Integral sich iiber die gesamte Ebene «a, f erstreckt, 
q(é, 7) = JQ, B) cosx € cos B y dx dB . (4) 


Substitution dieser Ausdriicke in die Gleichung (2) liefert mit der Abkiirzung 


N(a, B) = «4 D, + 20? B H + Bt D, (5) 
das Resultat : 


W(a, B) cosx x cos B y = [Q(a, 6) —- W(x, B) N(a, PY SKS) cosx& cosh y dé dy. (6) 


Im Sonderfall einer Auflast von der Form p = cos « x cos fy lieBe sich zeigen, daB die zugehdérige 
Bodeneinsenkung in der Form 


w = C(x, B) cosa x cos By (7) 
mit 
C(x, B) = [K (yx? + y?) cos x x cos B y dx dy (8) 
darstellbar ist?. Andererseits ist im selben Falle 
w = {K(s) cosa & cos B y dé dy. (9) 
E 


Aus dem Vergleich von (7) mit (9) folgt die Identitat 
poe) cos x € cos B y dé dy = C(x, 8) cosaxcospy. (10) 


1 Siehe z. B. K. Girkmann, Flachentragwerke, S. 300, Wien 1956, worin K , Ky an Stelle D,,, Dy steh 
2 Vgl. B.G. Korenev, Doklady Akad. Nauk, 78 (1951) §, 417, ao eames 


Aa ot ie Lt Cael oui 


bs sa AY CTT 


we 


pany ee ae ee 


Been S 


Pree 
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Setzt man nun den linken Teil von (10) in die Gleichung (6) ein, so lautet das Ergebnis 
W(x, B) = [Q(x, B) — N(o, 8) W(x, B)] Cla, B) « (11) 
Die Funktion W(q, 8) 148t sich hieraus bestimmen. Der Ansatz (3) liefert somit endgiiltig 


Q(a, B) C(a, B) 


w(x, y) = T+ C(a, A) Na, B) cosa xcosfh ydx dp, (12) 
E 
worin nach dem Fourierschen Theorem 
1 
O(a, B) = ree | q(, n) cos x € cos B n dé dy (13) 
E 


anzunehmen ist. Die zugehérige Bodenpressung ergibt sich schlieBlich aus der Gleichung (1) zu 


Pe.) = | O(c, B) cosa xcosPydadp. (14) 
E 


1 + C(a, B) N(@, B) 


Fiir den wichtigen Fall einer in € = 0, ” = 0 angreifenden Punktlast P gilt Q(a, 6) = P/4x?, 
und es ergibt sich dementsprechend eine elastische Flache der Platte von der Form 


Pe, Cas8) 
w ae | | ry cee by Nw H OO8 #008 BY ds a, (15) 
6 0 
wahrend sich die Bodenpressung nach 


vse) cos ax cos Py da dp (16) 
O70 


1+ C(@, B) N(@, B) 


verteilt. 
Von der Beziechung (8) schlieBt man noch allgemein auf die Reziprozitatsformel 
K(x? +?) = ae : C(x, B) cosa x cos B y dx dp . (17) 
s ; 
Mit 
raPtR, y= EP (1) 


14Bt sich die Relation (8) auch einfacher und zwar in der Form 
NY Fe) ENG) Sy) Ay (19) 
schreiben’, sofern J, die Besselsche Funktion nullter Ordnung bezeichnet. 


2. Sonderergebnisse in Abhangigkeit von den Eigenschaften der Unterlage. Je nach der Hypo- 
these, die zur Erfassung der mechanischen Eigenschaften der Plattenunterlage gewahlt wird, — 
kommt man im Einzelnen zu den folgenden Ergebnissen: 


a) Bodenziffer-Theorie. Der Winkler-Zimmermannschen Annahme gem4f ist C(y) eine 
von y unabhangige Konstante. Aus dem Vergleich der Ausdriicke (15) und (16) ergibt sich dann 


Cy) = = (20) 


wo die ,,Bodenziffer‘* k das konstante Verhaltnis zwischen dem jeweiligen Bodendruck und der 
betreffenden Einsenkung bedeutet. 


b) Elastischer Halbraum. Sind E, und », der Elastizitatsmodul und die Querdehnungszahl 
des isotropen Halbraumes, auf dem die Platte aufruhen soll, so gilt fiir den Kern K der Ausdruck 


Kr) =—,"8 (21) 


a Egor’ 


1 Wegen der Umformung 4hnlicher, auch im Abschnitt 4 vorkommender Integrale siche G. N. Watson, 
Theory of Bessel Functions, S. 124, Cambridge 1948. 
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und die Formel (19) liefert 


2(1— 7% 
CO) = aa : ue (22) 


c) Allgemeiner Fall. Setzt man allgemeiner 


B 
Kk) == (23) 
(a 
wo B und m wiederum Konstanten sind, so ergibt sich ganz entsprechend 
CC ae (24) 
d) Exponentialkern. Der Kern 
SS 
K(r) = 2 m2 é aim (25) 
mit den Konstanten B und m fihrt schlieBlich zu 
Cy) = 240 Begs (26) 


3. Sonderlosung fiir H = /D, Dy» Zu allererst soll dieser zu bedeutenden Vereinfachungen fih- 
rende Fall betrachtet werden. Bei D, =- D, sind die Steifigkeiten des Tragwerkes in den beiden 
Richtungen zwar verschieden groB, aber seine Torsionssteifigkeit entspricht etwa derjenigen eine 
rippenlosen Platte. Ursache der Anisotropie kann in diesem Falle eine in den Richtungen x und y 
verschieden starke Plattenbewehrung sein. In diesem Abschnitt soll nur die Bettungsziffertheorie 
beriicksichtigt werden, da nur sie, mathematisch gesehen, eine wesentliche Vereinfachung des 
Problems mit sich bringt. 

Mit den Bezeichnungen 


D=/D,D,, — 1=\/P (27) 
hat man zunachst 
' Doan D2. Bea. D,=D?. (28) 
Ks sei weiterhin 
xyA=%,, ee ae py. =6,, (29) 
und schlieBlich 
v1 =Voi + Bi. (30) 
Aus den Gleichungen (5), (28) und (30) folgt alsdann 
N(a, B) = Do? A> + PAP = Dy. (31) 


Kine Einzellast P greife nun in « = y = 0 an. Mit C(y) = 1/k geht dann der allgemeinere Aus- 
druck (15) fiir die zugehérige elastische Flache tiber in 


co © 


ets cos &, x, cos f, y 
eh Eat don dB, « (32) 


Dieser Ausdruck unterscheidet sich durch nichts von dem entsprechenden, nur in den Koordi- 


naten x, y, geschriebenen Ausdruck fiir die Durchbiegung einer schwimmenden isotropen Platte!. 
Setzt man noch 


Tp: ie | | 
l=\/>. Vey = 6s (33) 
so laBt sich (32) in folgender Weise durch die Thomsonsche Funktion kei ausdriicken?: 
ee OP Peat PP. ./x2 2 oie 
w=— zap kee — gap kell BA +9) (34) 


Vel. S. Woinowsky-Krieger, Ing.-Arch. 3 (1932) S. 246 und Ing.-Arch. 17 (1949) S. 142. 


2 Es gilt kei @ = (2/2) Im Y,(o '7!*) — (a/2) Re J,(o e'7/*) = — (2) Z.(o), wo J, ¥, B 
: : g K ny, : 
tionen und Z, eine der von F’, Schleicher tabulierten Le seis 6 nee Le 
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Héhenlinien dieser elastischen Flache sind in der Richtung der gréSeren Plattensteifigkeit lang- 


gezogene Ellipsen. Der Biegepfeil unter der Last ergibt sich mit kei 0 =— 7/4 fir 0 = 0 zu 
EF. 
(0) =—___.. 
ay 8 Vk D (35) 


Da Biegemomente einer diinnen Platte unter einer Punktlast stets unendlich werden, soll die 


Kraft P gleichmafig auf eine kleine Kreisflache vom Halbmesser c <1 mit dem Mittelpunkt in 


* = y = 0 verteilt werden. Die Funktion kei 9 mége nur mit den folgenden Anfangsgliedern ihrer 
Entwicklung angesetzt werden: 


: Qe a 2 
kei o ete 7 ot 9,209 0 : (26) 


Die Kriimmung 0?w/0x? der Flache (34) berechnet sich hiernach zu 


aw PA i ( 14 “ i) + ee oe bale = 0,616]. (37) 


Ox? 42D 


Indem man noch x = r cos 3, y = rsin # setzt, erhalt man im Mittelpunkt der angenommenen 


Druckflache 


wt 
0?w 


c 23 
Pi 1 : 
ae | r ar | Im = + = In (A cos? B + A+ sin? #) + 2? (A + A-1 tg? 8)? + 0,616 dj. 
é 


Nach Ausfithrung der Integration erhalt man die erste der beiden folgenden Formeln, und wenn 
man in ihr J durch 1/A ersetzt auch die zweite: 

0?w Pi l 2A (A— 1)? 
Ox? reo Ua eeeree cog 2 


Sh PALE dc 8 Od @— 1p 
te inp! (; aa +1 0,610. 


0,610 


Im Ansatz fiir die Biegemomente 


0?w 


0?w 0?w Ow 
me =—D. (Fa +% Be) my = —D, (52 +% 53) 


mu die Relation D,v, = D,¥,, die Querdehnungszahlen y, und y, betreffend, beriicksichtigt 
werden. Fiihrt man dementsprechend die Bezeichnung y, =v, v, =v/? mit v? = 7, v, ein, 
so gelangt man zu dem folgenden Ausdriicken fiir die Biegemomente im Zentrum der Druckflache: 


P(+») 2A aly 

m,(0) => eel In ( 2 22 as r) 2 aE 1 =f 0,616 ? (39) 
PA(1 +7) fae ® (A—1) 

m,(0) ==" n(, ergyee aoe + 0,616]. 


4. Struktur einer allgemeineren Lésung. Die Last sei zunichst iiber einen Kreisumfang vom 
Halbmesser c mit dem Mittelpunkt im Koordinatenanfang gleichformig mit der Intensitat p= P/27¢ 
verteilt. Unter Einfiihrung der Polarkoordinaten 


* 


E r=jeqy, ~ tea=X, 
|? + | - me 
ya=le+R, tep=“. 
ergibt sich nun auf Grund der Gleichung (13) das Resultat 
22 
Q(y, ¢) = £5 | cos (y ¢ cos g) cos (yc sin y) dp = f- Joly ¢) - (41) 


0 
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Der Ausdruck (5) 1aBt sich dementsprechend umschreiben in 
Ny, 9) = ® (42) 
mit 


@ = D, costy + 2 H sin? y cos? y + D, sinty . (43) | 


Die Gleichung (12) liefert somit die elastische Flache 


cos (y x cos @) cos (y y sin 7) | 
Soe J Sake. wi 1+ Cy) ® do . (44) 
0 
Ks sei nun weiter 
1 co 
ei 45 
1+ Cy) d— 2 Asal) cos2ng. ( ’ 


Bestimmt man die beiden ersten Koeffizienten dieser Fourier-Entwicklung durch irgendein Ver- 
fahren, so ergeben sich alle folgenden mit Hilfe einfacher Rekursionsformeln. Geht man an- 
schlieBend mit der Reihe (45) in den Ausdruck (44) ein, so hat man vorerst das Integral 


2a 


f cos (y x cos g) cos (yy sin g) cos 2n gy dp = 22 (—1)" Jenly 7) cos 2nd (46) 
0 
zu beachten. Man gelangt auf diese Weise zu dem Ergebnis 
300 0) a eA Ate LAY) OO) Jal) Jane (47) 


wobei C(y) nach wie vor nur von den Eigenschaften der Unterlage abhangig und durch den Aus- 
druck (19) gegeben ist. Wie der Vergleich von (12) und (14) lehrt, erhalt man die der Einsenkung 
(47) entsprechende Verteilung des Bodendruckes, indem man in (47) an Stelle von C(y) einfach eine 
Eins setzt. Im Sonderfall einer Punktlast P im Koordinatenanfang hat man nur p = P/2ac¢ 
anzunehmen und zum Grenzfall\c > 0 tiberzugehen. Das Resultat ist dann 


w(r, 9) = 5h SY (—1)t 008 2 0B f Analy) CC) Janly 17 ay (48) 


Im Fall D, = Dy also gleichberechtigter Richtungen x und y verschwinden, wie man sich leicht 
iiberzeugt, alle Glieder mit n = 1,3,..., und die vorstehenden Summen behalten nur Glieder 
mit Ay, A,, A, usw. 


5. Biegepfeil und Bodenpressung unter einer Punktlast. Wir beschranken uns hier auf die Angabe 
dieser zwei GréBen in einigen einfacheren Sonderfallen. Der allgemeine Ausdruck fir den Biege- 
ee wenn man in (15) x = y = 0 setzt und von den Bezeichnungen (40), (42) Gebrauch 
macht, 


x/2 fo) 

ay C(y) y dy 

oy aa we eet ay 
0 0 


Die maximale Bodenpressung erhalt man aus (49), indem man vorerst im Zahler des Integran- 
den C(y) durch eine Eins ersetzt und alsdann die Integration ausfihrt. 


a) Bodenziffertheorie. Fir C(y) = 1/k liefert die vorstehende Formel den Ausdruck 


7/2 fo) 1 


/2 
Ye d P d Fj 
w(0) saat dp | ——— = £ (50) 
I k+y*@ : 
A F Y Ax Vk F \@ 


Bei einem in den Richtungen xund y gleich steifem Balkenrost ist mit D, = D, = D, und H=7D, 
zu rechnen, worin 7 eine gegebene Zahl ist. Es ergibt sich hierfiir, der Gleichung (43) entsprechend, | 


D = D, (1 — é* sin? 2 ?) (51) ot 
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mit 
ef = — (52) 


und weiterhin, mit 2m = y und l, = /D,[k ; 


2 
w(0) a, = 0,125 a, —+. (53) 


% 
‘ox Iz dy a 
8x)kD, | Vl—esin?y 8YkD, D, 
0 
Hierin ist 
2 
0 =— K(e) (54) 


mit dem Modul des elliptischen Integrals ¢. Die Beiwerte o, sind in der Tabelle 1 zusammengestellt. 


Der maximale Bodendruck ergibt sich durch Multiplikation des Biegepfeiles (53) mit der Boden- 
ziffer k. 


b) Elastischer Halbraum als Unterlage. Mit den im Abschnitt 2 eingefiihrten Bezeich- 
nungen sowie den Abkiirzungen 


neces Gh D. Le 
ee ee e 
ergibt sich der Wert von C(y) zu 1/(y ko) und weiterhin, auf Grund von (49), der Biegepfeil gleich 


7/2 foe) n/2 


P dy 2 Py3 dp 
00) <a | dp | ee = Soe | ae: (9) 


0 0 0 


Indem man den vorstehenden Integranden durch C(y) = 1/(y ky) dividiert, erhalt man fiir die 
Bodenpressung unter der Einzellast den Ausdruck 


n/2 co 


7/2 
Pk y dy 2 Pk 73 [ dp 
OO bee eg | om or 
0 0 0 


Zwei Falle mégen nun speziell behandelt werden. 


1. Balkenrost mit D, = D, = D,, H =7 D,. Die Funktion @ ist in diesem Falle von der 
Form (51). Mit 2m = y erhalt man also , 


70/2 


2P 3 dy Py PI = 
= = - = 0,192 
(0) on hale D4} i Vie sint p 9 k2ls D8 Xs 0, 9 Xe D, ? ( ) 
6 


sofern ¢? = (1 — 7)/2 ist und a, den mit 2 multiplizierten Wert des vorstehenden Integrals be- 
zeichnet. Bei kleinen Werten ¢? berechnet sich dieses und ahnliche Integrale wohl am einfachsten, 


-indem man den Integranden nach Potenzen von ¢? sin? y entwickelt und die Reihe gliedweise aus- 
_ integriert. Bei gréBeren Werten ¢? kommt vor allem die Anwendung der Simpsonschen Formel 
in Betracht. 


Die Werte «, sind in der Tabelle 1 angegeben. Die durch den Verlust der Torsionssteifigkeit 
bedingte VergréBerung des Biegepfeiles erweist sich nach diesen Ergebnissen als weniger ausgepragt 
als bei der schwimmenden Platte. 

Die GréBe der maximalen Bodenpressung berechnet sich nach der Gleichung (57) zu 


n/2 = 
2 Pk2/3 13 dy i 21 is 59 
P(0) = aie | (—é sin? ys ( 4] gees a ©9) 
0 


9 x D3 aD 


-wobei sich x, ganz ahnlich ergibt, wie der Beiwert a, im Ausdruck (58). Auch die Zahlenwerte 
_ fiir xg sind in der Tabelle 1 zu finden. 


PPP eet ne ee ee EN 


2. Stahlbetonplatte mit D, = Dj, D, = D#?, H = D. Aus dem Vergleich von (31) mit (42) 
erhalt man in diesem Falle die Funktion 


® = D(A cos? y + Asin? y)? = DA [1 — (1— 7’) sin’ g}’ , 
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und die Formel (56) liefert 
n/2 oy 
2 PAzls V3 dp ea x PI 60 
w(0) = grep | [I—d—2)sn® pp? — One DIUS % Ne ear 2 (ey 
0 


Der Beiwert «, bezeichnet hier den mit 2 /?/5/ multiplizierten Wert des betreffenden Integrals, 
der wiederum identisch ist mit dem Integral der Gleichung (59). 3 
Die zugehérige gréBte Bodenpressung erhalt man ganz ahnlich aus der allgemeinen Formel (57). 


Es wird 


2/2 = 
2 P k3/8 als Y3- dp PS see ee iB 
P00) =~ DIB Toad aaiph® 9 (ee ed) 
0 


falls der Beiwert x, den mit 2 A4//z, vervielfachten Wert des Integrals bedeutet. Auch die Werte x, 
und «; finden sich in der Tabelle 1. 


Tabelle 1. Beiwerte «; zur Berechnung des Biegepfeiles und der Bodenpressung 


Fall: D, =Dy=D, H=nD, Fall: D, Dy = H?, 4 = (Dy/D,)1/4 
Bodenziffer- Elastischer Halbraum Elastischer Halbraum 
" Theorie a 
Oy Xe Xs OX, Os 
0,0 1,180 1,115 1,202 — — = 
0,1 1,155 1,099 1,215 0,1 0,753 1,442 
0.2 1,132 1,085 1,181 0,2 0,870 1,253 
0,3 1,110 1,072 1,151 0,3 0,925 1,148 
0,4 1,091 1,060 Lali) 0,4 0,955 1,092 
0,5 1,073 1,048 1,100 0,5 0,974 1,053 
0,6 1,057 1,037 1,076 0,6 0,986 1,029 
0,7 1,041 O27 1,055 0,7 0,993 1,014 
0,8 1,027 1,018 1,036 0,8 0,997 1,006 
0,9 1,013 1,009 1,018 0,9 0,999 1,002 
30 1,000 1,000 1,000 1,0 1,000 1,000 


Anmerkung: 0,(1/A) =a,(A);_ (1/4) = «;5(A) 


6. Biegemomente unter der Einzellast. Zur Berechnung der Biegemomente soll die Einzellast P 
von Anfang an auf eine Kreisflache von einem Halbmesser c verteilt werden. Liegt das Kreis- 
zentrum im Koordinatenanfang, so liefert die in Polarkoordinaten r, # statt in den Koordinaten é, n 
geschriebene Formel (13) 


c 2a 
00-9) =a aa | rar { cos (xr cos #) cos (fr sin #) dP 
6 6 


c 


ie) P 
=a [2 bly) rd =r hy) (62) 
6 
Die weitere Untersuchung mége auf den Fall eines ,,schwimmenden“ Balkenrostes beschrankt 
werden. Nimmt man nach wie vor D, = D, = Dy an, so geniigt es, zur Angabe der Biegemomente 
in der Mitte des Lastkreises den Wert 
: +00 + co 

ow ow Q(a, B) Cla, B) y? 

Aze(0) == as RE Oe ’ 2P)Y 
w(0) = (Gr + He ae [Aaa eS) 


—c —oc 


der sich aus (12) ableiten 1a8t, zu wissen. Geht man wiederum zu Polarkoordinaten iiber, so kommt 
mit Riicksicht auf (62) und auf die GréBen C = 1/k und N = yt © 


2% fee} 


P 2 Tye 
Aw(0) = —+—— peers: dy. (64) 
0 
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Setzt man weiterhin zur Abkiirzung 


ieee 
so ergibt sich 
22 
on P dp 
Aw(0) = = d et (66) 


0 


= | a | YP Sly ody [ Y? Sly 6) dy 
4 Q- —- Se Wak oo . 
pe apt : P+ (uy—i?? : Pt (wyi? 


Nun gilt fir —1 < Rv) < 2m+ 2 » K(k) > 0 allgemein! 


CO 
a TDS, (ae) ei Ko naib) 
(x? + 2)" 1 2™ m! j 
0 


In unserem Falle ist anzunehmen 
m=O, Ley bee a a= ¢, k=py—i bzw. k=pyi 
mit dem Ergebnis 


P= 5 4[—-Vi K (cui) +7 K (cus). (67) 


Die Macdonaldschen Funktionen K, (---) lassen sich ihrerseits durch die Thomsonschen Funktionen 
ersetzen®. Man erhalt nach einigen Vereinfachungen 


1 dkei(cp) _ 


1 y 
js A ee kei'(c u). (68) 


_ Nun ist auBerdem 


kei'(c us) il 
ex Re Fas In (cu) + 0,308 . 


Die Formel (66) liefert somit 
27 
P d 
Aw(0) = — 7 i [—In (cu) + 0,616] 2. (69) 
: 6 


Fir den Balkenrost gilt hier nach Gleichung (51) 
® = D, (1— é* sin® 2). 


Man kommt also mit der bereits friher benutzten Bezeichnung |, = 1D, /k zu einem Sonderwert fiir yu 


4 
= \s — 7 — ¢2 sin? 2 ele 
: Der Ausdruck (69) laBt sich hiernach wie folgt umschreiben: 


7/2 
12 I il j d 
2 Aw(0) => 72 D, | fn 5 + 0,616 + af In (1 — 6 sin? w) ay i (70) 
d 6 


Nun ist in Anbetracht von é? = (1 as y)/2 


7/2 
dy I 


1 — é? sin? p = y2a +7) Y 


0 


| 


Es gilt also im Endergebnis 
is Oe : h 4 0,616) + + 1| 
Aw(0) => xD, Faves; (In i Ie 0, ae 4 ’ 


£ 1G, N. Watson, Theory of Bessel Functions, S. 425. 
-: Siehe z. B. A. Erdélyi, Tables of Integral Transforms, Bd. 1, S. 371, New York 1954. 


\ 
Fe 


} ese 
4 ; 
eet. > 
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worin : 
7/2 
an In (1 — ¢? sin? y) 
=f 1—eé? sin’ yp ae 
6 
sein soll. 
Die Biegemomente m,(0) = m,(0) = — [D, (1 +») Aw(0)]/2 driicken sich somit aus in der 
Form 
m{0) = my(0) =~ (aint +B), (71) 
sofern man von den Bezeichnungen 
Tee 1 
«=| po und fp =0,616« +> 


Gebrauch macht. Die Beiwerte « und f sind in der Tabelle 2 zusammengestellt. 


Tabelle 2 

Beiwerte « und B zur Berechnung der Biegemomente 
n a B 1 & B 
0,0 1,414 0,453 0,5 L155 0,576 
0,1 1,348 0,490 0,6 1,118 0,588 
0,2 1,291 0,520 0,7 1,085 0,598 
0,3 1,240 0,543 0,8 1,054 0,606 
0,4 1,195 0,561 0,9 1,026 0,612 
0,5 DLS ois| 105000 1,0 1,000 | 0,616 


Auf die Berechnung der Biegemomente im Falle des elastischen Halbraumes als Plattenunterlage 
gehen wir nicht weiter ein. Es mége indes die folgende diesbeziigliche Bemerkung gemacht werden. 
Hat man es mit einer isotropen Platte zu tun, so erweisen sich die zugehérigen Ausdriicke fiir die 
Biegemomente unter einer stark konzentrierten Einzellast in beiden Fallen, das hei®t also bei der 
Anwendung der Bodenziffertheorie und beim elastischen Halbraum als Unterlage, als formal 
identisch!. Nur ist als Bezugslange im ersten Fall 


4 
D 
1=//2 
3 


D 
= |/ _— 
ko 
mit ky = E,/(2 — 2 »?) einzufihren. 

Héchstwahrscheinlich lieBen sich auch die Formeln (39) und (71) auf den Fall des elastischen 
Halbraumes mit einer fiir technische Zwecke ausreichender Genauigkeit anwenden. Nur ware, 
falls man so vorgeht, in der Gleichung (39) mit J, = /Dik, statt mit ] = | Dik und in der Glei- 
chung (71) mit J), = VD, /o statt mit 1, = /D,/k zu rechnen. 


und im zweiten 


(Eingegangen am 14. April 1959.) 


Anschrift des Verfassers: Professor Dr. S. Woinowsky-Krieger, Quebec (Canada), 
Université Laval Faculté des Sciences Boulevard |’Entente. 


1 S. Woinowsky-Krieger, Ing.-Arch. 17 (1949) S. 146. 
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Strémungsverhiltnisse in einem Diffusor 
mit vorgeschaltetem Kriimmer* 


Von K. H. Pohl 


1. Einleitung. Die vorliegende Arbeit soll zur Kenntnis der noch weitgehend unbekannten 
Strémungsverhiltnisse in einem Diffusor mit vorgeschaltetem Kriimmer beitragen. Da eine Lisung 
dieses Problems durch Integration der hydrodynamischen Bewegungsgleichungen bei dem augen- 
blicklichen Stand der Theorie aussichtslos ist, steht das Experiment ausgesprochen im Vordergrund. 
Als Hauptaufgabe ging es dabei um die Feststellung der auftretenden Verluste und die Aufnahme 
von Geschwindigkeitsprofilen. 

Um den Umfang der Versuche zu begrenzen, wurden fiir alle Messungen die Diffusor- und Kriim- 
mermaBe bis auf den Umlenkungswinkel des Kriimmers beibehalten. So konnte also nur eine 
grundsatzliche und keine umfassende Klarung der Strémungsvorginge erreicht werden; denn dazu 
ware es nétig gewesen, auBer dem Umlenkungswinkel des Kriimmers noch andere Parameter, z. B. 
das Seitenverhaltnis des Kanalquerschnitts, den Erweiterungswinkel und das Querschnittsverhaltnis 
des Diffusors und vor allen Dingen das Verhaltnis des mittleren Krimmungsradius zum hydrau- 
lischen Radius des Kriimmers zu verandern. Alle Versuche wurden mit Luft ausgefiihrt. 


2. Die Versuchseinrichtung und die Durchfiihrung der Versuche. Der Versuchskanal von 
quadratischem Querschnitt ist an die Saugseite eines Radialgeblases angeschlossen, wie dies Abb. 1 
zeigt. Die Strémung wurde in den Querschnitten IJ, III, IV und VII ausfihrlich mit MeBsonden 
untersucht. In den Querschnitten I, V und VI wurden nur Wanddriicke gemessen. Die fiir diese 
Querschnitte geltenden GréBen sind durch die Indizes 1 bis 7 gekennzeichnet. 
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Ausgleichstrecke 
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oa : Abb. 1. Versuchsanlage (in der Kanalanordnung A III). 


Die Luft tritt durch eine Einlaufdiise mit ebenfalls quadratischem Querschnitt in den Kanal ein. 
Ein Gleichrichter und eine auf diesen folgende Kontraktion der Hinlaufstrecke im Verhaltnis 5:1 
_ sorgen fiir die Beseitigung der Einlaufstérungen. ee ¢ 

Die folgende gekriimmte Kanalstrecke besteht wahlweise aus einem oder zwei Kriimmerstiicken 
- mit einem Umlenkungswinkel 6 von je 90°, so daf mit diesen beiden Teilen eine Umlenkung um 90 

oder 180° méglich ist. Die Kriimmer haben eine Querschnittsflache von 0,18 < 0,18 m2; das Ver- 
hiltnis des mittleren Kriimmungsradius 7,, zum hydraulischen Radius ry ist r[lq = 9,09. 


* Gekiirzte Fassung der vom Verfasser im Jahre 1954 im Max-Planck-Institut fir Strémungsforschung 
| Gn Géttingen durchgefiihrten Diplomarbeit (Ref.: Prof. Dr. Albert Betz). 
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Der unmittelbar an den Kriimmer angeschlossene Diffusor hat einen halben Erweiterungs- 
winkel von 5° und ein Flachenverhaltnis von Ein- und Austrittsquerschnitt von F3/Fy = 0,4, also 
eine Lange von 0,6 m. 

Das dem Diffusor folgende Kanalstiick ist 2m lang und hat eine Querschnittsflache von 
0,285 X 0,285 m?. Diese ,,Ausgleichsstrecke* soll zur VergleichmaBigung des sich am Ende des 
Diffusors ausbildenden Geschwindigkeitsprofils dienen und so eine teilweise Umsetzung von kine- 
tischer Energie in Druck erméglichen. 

Nach Verlassen dieser Ausgleichsstrecke stromt die Luft durch das von einem Drehstrommotor 
angetriebene Geblase. Der Durchflu8 konnte mit einer Drosselklappe am Ventilator geregelt wer- 
den. Bis auf die Einlaufstrecke, die ganz in Blech ausgefiihrt ist, bestehen die Seitenwande der 
quadratischen Versuchsstrecke aus Blech, Boden und Deckel aus Sperrholz. Der Diffusor hat 
hélzerne Seitenwainde, Deckel und Boden sind aus Glas. So kann man die Strémung im Diffusor 
qualitativ mittels Faden oder Rauch verfolgen. 

Durch Veranderung des Umlenkungswinkels 6 des Kriimmers ergeben sich die folgenden drei 
Kanalanordnungen A, an denen die Messungen vorgenommen wurden: 


AI ;:6= _ 0°: der gerade Kanal mit Diffusor, 
AII :6= 90°: Kanal mit einem 90°-Kriimmer vor dem Diffusor, 
A III: 6 = 180°: Kanal mit einem 180°-Kriimmer vor dem Diffusor. 


Die Durchfiihrung der Versuche erforderte die Messung des Staudruckes q = 9/2 v?, des 
Gesamtdruckes g und des statistischen Druckes p in den Querschnitten IJ, III, IV und VII, 
dabei konnte natiirlich je einer der drei Driicke aus den beiden anderen nach der Definitionsglei- 
chung p + q =g berechnet werden. AuBerdem wurden an simtlichen Mefstellen I bis VII die 
statischen Driicke an der Wand gemessen. In allen drei Kanalanordnungen AI, A II und A III 
wurde die Strémung in Abhangigkeit von drei Re-Zahlen untersucht. 

Die Messungen wurden mit einem unter 0° gegen die Kanalachse eingestellten Staugerat aus- 
gefiihrt. Das ist méglich unter der Annahme, dab die Strémungsrichtung nur wenig von der Kanal- ~ 
achse abweicht, so daf man, ohne einen merklichen Fehler zu begehen, den Absolutwert v der 
Strémungsgeschwindigkeit, den das Staurohr bis zu Abweichungen von 15° richtig anzeigt, fiir die 
in die Axial-Richtung fallende Komponente v, setzen kann. An einigen Punkten des Kanals wurde 
diese Annahme mit einer Vier-Loch-Sonde nachgepriift, und, solange die Strémung nicht abgerissen 
war, bestatigt gefunden. Trat irgendwo Riickstrémung auf, so wurde der Staudruck mit dem 
um 180° gedrehten Staugerat ermittelt. 


3. Die Versuchsergebnisse. Alle MeBergebnisse und die aus ihnen abgeleiteten Werte wurden 
mit einer geeigneten GréBe dimensionslos gemacht. Die Darstellung hat den Vorteil, daB die so 
erhaltenen Ergebnisse fiir geometrisch ahnliche Kanale und beliebige inkompressible Flissigkeiten 
gelten, wenn man nur darauf achtet, daB in den beiden zu vergleichenden Fallen die Re-Zahlen 
ungefahr gleich gro sind. In diesem Fall definiert man die Re-Zahl zweckmafig durch 


Re 40H | (1) 


mit y = kinematischer Zahigkeit, v, = mittlerer Geschwindigkeit im engsten Querschnitt IJ und 
To = 4,/4 = hydraulischer Radius des Querschnitts IJ. Im folgenden werden die dimensions- 
losen GréBen im Unterschied zu den dimensionsbehafteten mit einem Strich gekennzeichnet. 
a) Die Geschwindigkeitsverteilung. Die Geschwindigkeiten wurden jeweils mit der mitt- 
leren Geschwindigkeit 1, des betreffenden Querschnitts dimensionslos gemacht: 
ey Ne v,(x' y') at 
| AS ae opersermis (2). 
In den Querschnitten IJI und IV traten zum Teil zeitlich stark schwankende Driicke auf. An 
solchen Stellen liefert die Messung der Strémungsgeschwindigkeit mit dem Staugerdt unvermeidlich 
zu hohe Geschwindigkeitswerte ; denn man mift mit dem Staugerat iiber eine genugend lange Zeit T 


den Mittelwert des Geschwindigkeitsquadrates v2. Bekanntlich ist aber, wenn man v = 0 + v(t) 
setzt, 


re Te eee 3 
faU tg [Ha> o. (3) | 
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pp ain Querschnitt II v(t) a 0 ist, ist die dort aus v2 ermittelte mittlere Geschwindigkeit v, der- 
ee Wert, welcher der wirklichen mittleren DurchfluBgeschwindigkeit am nachsten kommt. Mit 
resem ergeben sich aber die mittleren DurchfluBgeschwindigkeiten in den anderen Querschnitten 
auf Grund des iiberall gleichen Volumenstromes V. Um die langwierige Ermittlung von 3, bei 
. td 
jedem Versuch zu ersparen, wurde der Volumenstrom V durch Messung des Wanddruckes p, un- 
mittelbar hinter dem Einlauftrichter unter Beriicksichtigung des Diisenkoeffizienten bestimmt. 
Die Verteilung der Geschwindigkeiten in den einzelnen MeBquerschnitten — dargestellt durch 
die Isotachen — bei den Strémungen durch die Kanale AI, AIT und A III und der Re-Zahl 5 - 105 


Abb, 2. MeBstelle ITT. Abb. 3. MeBstelle IV. 


ist in den Abb. 2 bis 12 wiedergegeben. Die Er- 
gebnisse bei anderen Re-Zahlen unterscheiden 
sich nicht stark von den dargestellten und zeigen 
daher nichts grundsatzlich anderes. 


: a ee : 58 
Alle Querschnitte sind in Strémungsrichtung ~ 


gesehen. Die Kriimmerinnenwand liegt links, die 
AuBenwand rechts. An Stelle der wirklichen 
Koordinaten x und y sind ihre mit der jeweiligen 
Kanalbreite (= Kanalhéhe) a, dimensionslos ge- 
machten Werte x’ = x/a, und y’ = y/a, verwandt. 
Die Querschnitte erscheinen dadurch stets in 
gleicher GréBe. Die Zahlen an den Isotachen be- 
deuten die dimensionslosen Geschwindigkeiten 
v(x’, y’) nach (2). Die Erstreckung der Kanal- 2 
achse vom Eintrittsquerschnitt IJ in den Kriim- 
mer in Strémungsrichtung ist als z-Achse ge- 
wahlt. Sie ist mit der Kanalbreite im engsten , 7 ai an Po ; 
Querschnitt II dimensionslos gemacht: z’ =2z/a,. 
An Hand dieser Abbildungen kann man folgen- 
_des erkennen: 
1. Kanalanordnung A I. Die Abb. 2 bis 4 zeigen die schon seit langem grundsatzlich bekannten 
Formen der Geschwindigkeitsverteilung im geraden Kanal mit Diffusor. In Abb. 2 ist deutlich 
die in jedem geraden Kanal von nicht kreisférmigem Querschnitt bei turbulenter Strémung auf- 
tretende Sekundarstrémung zweiter Art, die eine Ausbuchtung der Isotachen in die Kanalecken 
bewirkt, zu erkennen. 

Abb. 3 zeigt folgendes: Eine Betrachtung der Bernoullischen Gleichung lehrt, da die am Dif- 
fusor auftretende Verzégerung der Strémung eine Zuspitzung des Geschwindigkeitsprofiles nach 
der Kanalmitte zu bewirkt (vgl. Abb. 2 und 3). Wie bekannt, darf der Diffusor, wenn die Strémung 
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Abb. 4. MeBstelle VIT. 
Abb. 2 bis 4. Isotachen bei Kanalanordnung A 1. 
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nicht merklich abreiBen soll, im allgemeinen keinen gréBSeren Gesamterweiterungswinkel als 8° 
bis 10° haben. Der vorliegende Versuchsdiffusor mit einem Erweiterungswinkel von 10° versetzt 
die Strémung gerade in den ,,kritischen Ablésungszustand“. Schon eine geringe Stérung der Ge- 
schwindigkeitsverteilung kann Ablésung hervorrufen. Diese wird hier durch eine kaum merkliche 
Unsymmetrie in der Form des Diffusors verursacht. Die Strémung reiBt an der Oberseite des Dif- 


Abb. 5. Mefstelle IT. Abb. 6. MeBstelle III. 
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Abb. 7. MeBstelle IV. Abb. 8. MeBstelle VII. 


Abb. 5 bis 8. Isotachen bei Kanalanordnung A II. 


fusors ab, und es tritt Riickstrémung ein (schraffiertes Gebiet in Abb. 3). Die vorher sehr gute 
Symmetrie zur horizontalen Mittelebene y’ = 0,5 ist gestért. 

Am Ende der Ausgleichstrecke (Abb. 4) haben sich die starken Geschwindigkeitsunterschiede 
im Austrittsquerschnitt des Diffusors wieder einigermaBen ausgeglichen. Die Symmetrie des 
Profils zur Ebene y’ = 0,5 ist aber noch lange nicht wieder erreicht. 


2. Kanalanordnung A Il. Dagegen zeigen die Abb. 5 bis 12 fast ausnahmslos recht gute Sym- 
metrie zur Ebene y’ = 0,5 und gar keine Symmetrie mehr zur vertikalen Mittelebene x’ = 0,5. 
Das ist der Einflu8 des Kriimmers und der in ihm auftretenden Sekundarstrémung. Schon vor 
dem Eintritt in den Kriimmer (Abb. 5) wird die Strémung durch ihn in der Art beeinfluBt, daB sie 
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sich auf die Geschwindigkeitsverteilung, die fiir die Bewegung einer idealen Fliissigkeit im Kriimmer 
zu erwarten ist (vr = konst.), einstellt. Eine Verzégerung der Strémung an der rechten Seite 
des Kanals ist deutlich zu erkennen. Abb. 6 zeigt die Geschwindigkeitsverteilung im Austritts- 
querschnitt des Kriimmers, die durch Uberlagerung der Haupthewegung parallel zur Kanalachse 
und der Sekundarstrémung entstanden ist. Riickstr6mung an der Innenwandung des Kriimmers 
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Abb. 9. MeBstelle ITI. Abb, 10. MeBstelle IIT. 
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Abb. 11. MeBstelle IV. Abb. 12, MeBstelle VII. 
Abb. 9 bis 12. Isotachen bei Kanalanordnung A III. 


i i i ist die kinetische Energie der an der linken 
urde dort noch nicht gemessen. Erst im Diffusor ist die kinetise : g 
Wand flieBenden Fliissigkeitsteilchen so klein geworden, da die Strémung den dort bestehenden 
Druckanstieg nicht mehr iiberwinden kann und zuriickflieBt (Abb. 7). Das Totwassergebiet, das 
bei der Anordnung A I an der Oberseite des Diffusors entstanden war, kann sich hier unter dem 
Einflu®B der Sekundarstrémung nicht ausbilden. 
Abb. 8 zeigt dann, entsprechend wie bei der Anordnung A I, wieder einen Ausgleich der starken 
Geschwindigkeitsunterschiede im Vergleich zu Abb. 7. : 
Auf die Unterschiede der Kanalanordnung A IIT (Abb. 9 bis 12) gegeniiber der Anordnung A II 
soll im Abschnitt d) (Deutung der Versuchsergebnisse) eingegangen werden. 
3* 
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b) Der statische Druck p. Der Druckanstieg im Diffusor ist wegen der unvermeidlichen 
Verluste immer kleiner als er sich nach der Bernoullischen Gleichung fiir eine reibungslose Fliissig- 
keit ergibt. Ist der Erweiterungswinkel des Diffusors so groB, daB Ablésung der Strémung eintritt, | 
so liegt der Gewinn an statischem Druck weit unter dem nach Bernoulli errechneten. In der dem | 
Diffusor folgenden Ausgleichstrecke erfolgt aber durch eine VergleichmaBigung des Geschwindig- | 
keitsprofils ein zusatzlicher statischer Druckanstieg. So zeigt die Anordnung Diffusor—Ausgleich- | 
strecke wieder ein giinstigeres Bild. Ist p, der gemessene mittlere statische Druck im Querschnitt | 
k und p’;, = px/(o/2) 03 sein dimensionsloser Beiwert, dann ist Api, = pi’ — pr der gegeniiber dem | 
idealen nach Bernoulli berechneten statischen Druck p%’ bis zum Querschnitt k aufgetretene dimen- | 
sionslose Druckverlust. Demnach stellt die Differenz Api; = Ap, — Ap; den statischen Druck- | 
verlust gegeniiber dem nach der Bernoullischen Gleichung errechneten im Diffusor allein dar. 


Lei ie 


Kréimmer Difusor Ausgleichstrecke ] 


IV : VI 


Me&stelle 


Abb, 13. Verlauf der statischen Driicke. 


In den Abb. 13, 14 und 15 sind die Druckdifferenzen p,— p2, Ap; und Ap43 graphisch dar- 
gestellt. (Da der Versuchskanal nur vier MeBquerschnitte hat, war die Feststellung des genauen 
Kurvenverlaufes nicht méglich. Die in diesen Abbildungen gewahlte Darstellungsart der Kurven 
durch gerade Linienstiicke wird also dem wahren Kurvenverlauf nur annahernd entsprechen. 


c) De r Energiestrom. Ahnlich wie die Geschwindigkeitsverteilung kann man auch in jedem 
Querschnitt die Verteilung der durchflieBenden Gesamtenergie auftragen. Ist g = p + (0/2) v? 
der Gesamtdruck, so ist die durch ein Flachenelement df flieBende Gesamtenergie 


dE, = &, % af (4) 
und mit den dimensionslosen Werten gi, = g,/(0/2) v3, vj; = v,/0, und df’ = df/F, ihr dimensions- 
loser Wert 

, dE; 
Pas a ee (dob co? 
ivi Ge ae (5) 


Die sich hierbei ergebenden Kurven sehen natiirlich qualitativ ganz ahnlich aus wie die Isotachen. - 
Durch Integration tiber den ganzen Querschnitt erhalt man dann den durch diesen Querschnitt 
flieBenden Energiestrom. Er nimmt natiirlich in jeder Strecke um den in dieser Strecke auftreten- 


Be) 
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den Energieverlust ab. In Abb. 15 ist der Verlauf dieser Energieverluste ahnlich wie der der Druck- 


verluste in Abb. 14 dargestellt, und gibt eine gute Ubersicht iiber die Verluste in den einzelnen 
Strecken. 


d) Deutung der Versuchsergebnisse. Aus Abb. 15 ersieht man, daf die Verluste an Ge- 
samtenergie um so groBer werden, je starker-die Ablenkung durch den Kriimmer ist. Dies war ja 
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Abb, 14. Verlauf der Druckverluste. 
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Abb. 15. Verlauf der Energieverluste, 


auch zu erwarten. Aber aus einem Vergleich der Abb. 14 und 15 zeigt sich, daB die Verteilung der 


Gesamtenergie auf Druck- und kinetische Energie bei den verschiedenen Anordnungen doch recht 


 verschieden ist. Da andererseits die Verwertung der kinetischen Energie nur beschrankt méglich 


7A 


"ist, interessiert auch der Anteil der Druckenergie allein, da dieser auf alle Falle verwertbar ist. 


Hierbei fallt auf, da® der Druckverlust im Diffusor durch den davorliegenden Kriimmer nicht ver- 


- gréBert, sondern verringert wurde, insbesondere bei Anordnung A II. Man ersieht dies besonders 


38 K.H. Pohl: Stromungsverhaltnisse in einem Diffusor mit vorgeschaltetem Kriimmer Ingenieur-Archiy 


deutlich aus der Darstellung in Abb. 16, in der die Druckverluste Api3-gesondert dargestellt sind. 

Man kann diese Beobachtung durch folgende Uberlegung erklaren: 

Auf einen Fliissigkeitsstrom in einem gekriimmten Kanal wirken Zentrifugalkrafte ein, die 
eine Strémung senkrecht zur Kanalachse verursachen. Diese sogenannte Sekundarstrémung bildet 
sich als Doppelwirbel aus. Die resultierenden Stromlinien sind also zwei gegensinnig laufende 
Spiralen. Ist v die Geschwindigkeit parallel zur Kanalachse und r der Kriimmungsradius des 
Kriimmers an der betrachteten Stelle, so ist die auf ein Flissigkeitsteilchen wirkende Zentrifugal- 
kraft proportional v?/r. Daher ist es sicher berechtigt, v?/r als Ma® fir den Geschwindigkeitsbetrag 
oder die ,Starke‘ der Sekundarstrémung an dieser Stelle anzusehen. Nach der Potentialtheorie 
erhalt man die Geschwindigkeitsverteilung im Kriimmer v = const/r. Die aus der Kontinuitats- 
gleichung folgende mittlere Geschwindigkeit sei v. Bildet man nun den Mittelwert v?/r tiber einen 

Kanalquerschnitt fiir die beiden Geschwindigkeitsprofile v und », 

Ro=5-19° so ist v2/r > v2/r. 

APs — — Re=3,3-70° Das in einer Potentialstrémung an der Kriimmerinnenwand 
auftretende Geschwindigkeitsmaximum verschiebt sich in einer 
reibenden Strémung mit gréRer werdendem Umlenkungswinkel 6 
des Kriimmers in Richtung wachsender r. Der Betrag der Maximal- 
geschwindigkeit wird dabei gleichzeitig kleiner (vgl. Abb. 5 und 9). 
Andererseits kann sich kurz nach dem Eintritt der Strémung in 
den Kriimmer noch kein starkes Maximum der Geschwindigkeit 

| § ausbilden. Bei einem bestimmten Umlenkungswinkel 6, nimmt 

0 7) 780° och ; ; : ss c 

AI AL AZ also das Geschwindigkeitsmaximum seinen gréBten Wert an. Bei 

Kanalanordnung jedem anderen Winkel ist die Geschwindigkeitsverteilung gleich- 

férmiger. Da nun v2/r > v?/r ist, wird bei dem Winkel 6) die 

durch Integration von v?/r tiber den Endquerschnitt des Kriim- 
mers ermittelte Starke der Sekundarstrémung am gréBten. 

Beim Eintritt einer Kriimmerstrémung in einen Diffusor wird die kinetische Energie der Se- 
kundarstrémung, die einen wesentlichen Teil der Kriimmerverluste ausmacht, zum groBem Teil 
wieder in Druck umgesetzt. Da die Sekundarstrémung auBerdem die Totwasserbildung im Diffusor 
vermindert, werden auch die im Diffusor auftretenden Verluste kleiner. Da die Starke der Sekun- 
darstrémung bei dem Umlenkungswinkel 6, am gréBten ist, durchlaufen die Diffusorverluste ein 
Minimum. Die genaue Bestimmung dieses Winkels 5, war mit der vorhandenen Versuchsanordnung 
leider nicht méglich. 

Sind also in einer Kanalanordnung ein Kriimmer und ein Diffusor hintereinandergeschaltet, 
80 kann der Verlust in der gesamten Anlage wesentlich kleiner als die Summe der Einzelverluste 
sein. : 


A 


Abb. 16. Druckverlust im Diffusor AP 43 


4. Zusammenfassung. In der vorliegenden Arbeit wird die Strémung in einem Diffusor mit 
vorgeschaltetem Kriimmer untersucht. In dem Diffusor von quadratischem Querschnitt und einem 
Gesamterweiterungswinkel von 10° tritt bei einer Kanalanordnung ohne vorgeschalteten Kriimmer 
Ablésung ein. Die im Diffusor auftretenden Verluste sind dadurch betrachtlich. Das Ergebnis 
der Versuche zeigt, daB ein vorgeschalteter Kriimmer die Druckverluste im Diffusor vermindern 
kann. Dadurch wird der Druckverlust in der Gesamtanordnung kleiner als die Summe der Einzel- 
verluste von Kriimmer und Diffusor. Das ist auf den Einflu8 der sich im Kriimmer ausbildenden 
Sekundarstrémung zuriickzufiihren. Eine iberschlagige theoretische Betrachtung und eine genauere 
experimentelle Untersuchung zeigt weiter, daB die Druckverluste im Diffusor mit wachsendem — 
Umlenkungswinkel des Kriimmers ein Minimum durchlaufen. Das hangt damit zusammen, daB 
die Sekundarstrémung bei einem bestimmten Umlenkungswinkel am starksten ist. 


(Eingegangen am 23. April 1959.) 


Anschrift des Verfassers: Dipl. Phys. K. H. Pohl, Hamburgische Schiffbau-Versuchsanstalt, 
Hamburg 33, Bramfelderstr, 164. 
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Uber die Erfassung der Widerstandserhéhung im Seegang 
durch Energiebetrachtungen 


Von K. Eggers 


1. Einleitung. Vor kurzem erschien von Havelock! eine Studie iiber die Abhangigkeit der 
Dampfung der Stampfschwingungen eines schiffsaihnlichen Kérpers von Frequenz und Fahrt. 
Havelock untersucht eine Stromung, die er aus einer approximativen Quell-Senken-Verteilung 
erzeugt und ermittelt die Dampfung durch Integration des Drucks auf einer Flache, welche der 
mittleren Lage der Bodenflache des Kérpers entspricht. 

Fur die hier betrachtete Strémung laBt sich nun die Energie beechnen, die infolge der Wellen- 
bildung wahrend einer Periode abwandert?; es ergibt sich, daB ihr Wert der Arbeitsleistung der von 
Havelock angegebenen Dampfungskraft entspricht. Ebenso la@t sich aber auch die Ausstrahlung 
von Impuls in diesem Wellensystem berechnen; ihr Zeitmittel verschwindet nicht und kann somit 
nur aus einer Erhéhung des Wellenwiderstandes des stampfenden Schiffes in Fahrt herrihren, zu 
dessen Uberwindung ebenfalls Arbeit geleistet werden mu8. Bringt man diesen Anteil von der 
gesamten Energieabwanderung einer Periode in Abzug, so reicht dann der Restbetrag der Energie 
nicht mehr aus, um den Leistungen von Havelocks Dampfungskraften zu entsprechen; die Abwei- 
chung wachst mit zunehmender Fahrt und Schwingungszeit. 

Die sich hier offenbarende Diskrepanz zwischen der Druckintegrationsmethode und der Wellen- 
energiemethode verschwindet jedoch, wenn wir der gewahlten Quell-Senken-Verteilung zur Er- 
zeugung der Strémung eine Dipolverteilung von oszillierender Intensitat superponieren, welche die 
Variation des Deplacements fiir den Fahrtwiderstand erfaBt und damit die Approximation der 
Strémung auf der Kérperoberflache verbessert. Jetzt liefern beide Methoden tbereinstimmende 


_ Ergebnisse. 


aN os 


Die hier angewandte Verbesserung wurde fiir idealisierte Schiffsformen (Michell-Schiffe) bereits 
von Hanaoka’ eingefiihrt. Sie kann aber‘ allgemein fiir alle Kérper angesetzt werden, deren Poten- 
tial aus einer Quellbelegung proportional der Normalgeschwindigkeit der Strémung auf der Kérper- 
oberflache angenahert wird. Fir Ellipsoide in Schwingungen kleiner Amplitude gilt sie damit 
sogar streng*, so lange der Einflu8 der freien Oberflache auf die Strémung am Kérper vernachlassigt 
werden kann. 

Fiir die Energiemethode zur Bestimmung der Dampfungskrafte von Schwingungen in Fahrt 
wird die GréBe der Impulsabstrahlung bendtigt, aus der sich die Widerstandserhéhung infolge 
erzwungener Schwingungen ergibt; dariiberhinaus aber gestattet sie auch die Berechnung der 
Widerstandserhéhung in einem homogenen Wellensystem, das diese Schwingungen erzeugt. Dazu 
ist es keineswegs erforderlich, den Bewegungsverlauf explizit in seiner Phasenlage und in seinem 
Amplitudenverhaltnis zur ankommenden Welle zu kennen, soweit er nicht in die Quellverteilung 
zur Erzeugung des Potentials der Wellenreflexion eingeht; es geniigt die Kenntnis der Schwingungs- 
amplitude. Durch die Forderung, daf jetzt die gesamte abgestrahlte Energie dem gesamten ab- 
gestrahlten Impuls entsprechen mu, daB also die impulslose Leistung der Schwingungsdampfung 
vom entsprechenden Anteil des erregenden Wellensystems kompensiert wird, ergibt sich ein ein- 
facher Zusammenhang zwischen Dampfungsleistung und Widerstandserhéhung durch Fahrt gegen 
die erregende Welle, welche zur Widerstandserhéhung durch erzwungene Schwingung hinzugefiigt 
werden mu&. Die folgenden numerischen Rechnungen zeigen tibrigens, dafs beide Widerstands- 
anteile von gleicher Gré®enordnung sind, auch wenn die Wellenreflexion am festgehaltenen Schiff 
klein bleibt, und daB ihre Summe stets positiv ausfallt. 

Es ergibt sich damit die Méglichkeit, im Rahmen einer linearen Theorie die bisher nur einzeln 


- diskutierten Anteile der Widerstandserhéhung — Wellenreflexion, Widerstand in erzwungener 


Schwingung, drifting force — zusammen mit bisher noch nicht erfaften Termen gleicher GréBen- 


T. H. Havelock, Trans. Inst. Nav. Arch. 99 (1958) S. 131. London 1957, 

K. Eggers, Schiffstechnik 21 (1957) S. 140. ; ‘ 

T. S. Hanaoka, Proceedings of the Symposium on the Behaviour of Ships in a Seaway, Wageningen 1S 
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ordnung zu einem analytisch iibersichtlichen Ausdruck zusammen zu fassen. Bei vorgegebenen 
Amplituden der Schwingungen verschiedener Freiheitsgrade miissen hier nur noch die korrespon- 
dierenden Singularitatenverteilungen gefunden werden, was ein Problem der numerischen Analysis 
ist. Es ist uns damit vor allem méglich, durch Betrachtung charakteristischer vorgegebener Sin- 
gularitatenverteilungen die qualitative Abhangigkeit von Parametern der Fahrt und Frequenz 
zu studieren. Da®B der Zusammenhang zwischen Singularitéten und erzeugter Kérperform bzw. 
Korperbewegung dabei relativ offen bleibt, darf im augenblicklichen Stadium der theoretischen 
Entwicklung nicht ins Gewicht fallen. Zwar liegen schon Arbeiten von Haskind* vor, welche fiir 
Kérper von genau festgelegten Wasserlinien und Spantformen zahlenmafige Ergebnisse liefern — 
es fehlt hier jedoch im Singularitatenansatz der Dipolterm fiir den FahrteinfluB; die Widerstands- 
erhéhung der erzwungenen Schwingung, auf welche sich die Rechnung beschrankt, fallt deshalb 
negativ aus. Exaktere Rechnungen von Hanaoka und Maruo? (letztere unter Beriicksichtigung 
des Seegangs) beschranken sich auf einen sehr begrenzten Fahrt- und Frequenzbereich. Auch ist 
in diesen Rechnungen iibergangen worden, da8 Dampfung und Widerstandserhéhung iber alle 
Grenzen wachsen, wenn sich das Verhaltnis der Schiffsgeschwindigkeit zur Phasengeschwindigkeit 
einer Welle der betrachteten Frequenz dem Werte 1/4 nahert. 


Wahrend fiir erzwungene Schwingungen sich dieses Anwachsen schon aus den Versuchen von 
Golovato® erkennen lieB, zeigen neuere Experimente* in Ubereinstimmung mit unserer Analysis, 
daB der Effekt auch bei Fahrt eines freibeweglichen Modells gegen einen regelmaBigen Wellen- 
zug auftritt, und zwar nicht nur fiir die Froudesche Zahl % = 0,08, wo dieses Verhaltnis sich 
einstellt, falls Schiffslange und Wellenlange zusammenfallen, sondern auch bei % = 0,04, wo 
offensichtlich fiir die durch die Schwingung angeregte sekundare Oberschwingung das Geschwindig- 
keitsverhaltnis 1/, vorliegt. Aus den folgenden Berechnungen ergibt sich, da dieser Effekt auch 
bei Fahrt mit der Welle zu erwarten ist, und zwar einmal, wenn die Welle mit zweifacher Schiffs- 


geschwindigkeit voranschreitet, zum andern bei einer Schiffsgeschwindigkeit, welche gréBer als 
die der Welle ist. 


2. Geschwindigkeitspotential. Bei der analytischen Behandlung der Widerstandserhéhung im 
Seegang ist es iiblich, den schwingenden Schiffskérper zu ersetzen durch ein System von Quellen 
und Senken oszillierender Intensitat, welche nur die Translationsbewegung des Kérpers mitmachen. 
Die Ermittlung der Dampfungsenergie und Widerstandserhéhung stellen wir im folgenden vorerst 
nur dar fiir den Fall, da die periodische Stérung des Wellensystems bereits durch eine einzige 
Quelle dargestellt werden kann, solange der FahrteinfluB vernachlassigt wird; fiir kontinuierliche 
Quellverteilungen ergeben sich daraus wegen der Linearitat der Ansatze durch Integration analoge 
Ausdriicke, welche noch die Interferenzen zwischen verschiedenen Quellsystemen erfassen. 


Das Geschwindigkeitspotential y der Strémung relativ zu einer oszillierenden Quelle, welche sich 
mit Geschwindigkeit V parallel zur freien Obertlache in Tiefe f bewegt, wurde von Haskind, Brard', 
Hanaoka und Havelock auf verschiedenem Wege hergeleitet. Bei einer Quellergiebigkeit e‘”* lautet 
es im Aufpunkt r 


mm 
/ 


goa Lao cate arp + cdia are) 
R R Tt (K cos O— »/V)? T (Kens O- air) pe 


dG, 


eivt 


aan Ae, 


i Ky ett F sec? O athe 
ee KK (6 K, e %(' £9) + 6, K, oB(-Fo) | d0 (1) 


Dabei ist r der Vektor vom an der freien Oberfliche gespiegelten Quellpunkt zum Aufpunkt, 
R = |r| der Abstand vom gespiegelten Quellpunkt zum Aufpunkt, R der Abstand vom Quellpunkt 
zum Aufpunkt. Fiir ein rechtwinkliges Koordinatensystem (x, y, 2) wahlen wir die Richtungen: 
x = Fahrtrichtung voraus, y = backbord, z = nach oben. Ferner ist Eg ein Vektor mit Kom- | 


ponenten (i cos 0, isin@, 1), Eg der dazu konjugiert komplexe Vektor; © entspricht dem Aus- 
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breitungswinkel und K der Wellenzahl eines Systems von Wellen, iiber die integriert wird; im 
Integral der zweiten Zeile sollen K und @ miteinander verkniipft sein durch 


K, = K,(6) = K, cect @(! BE 1 Ee (1b) 


~ 2 2 
So Aal@) —K r 7 faa) ; oo 
dabei ist gesetzt: 


Ky = Le 
~ yp 
ee 
(ge ie Ve 
g K, 


SchlieBlich sind die Koeffizienten 6, und 6, so definiert, daB sie in dem @-Bereich verschwinden, 
in dem K, und Ky nicht reell ausfallen, im iibrigen gilt 


,=+1 fir [ol<<, | 


\ falls K,, Ky reell. 


§=—1 fir joj>%, | 
ee tir 0 = (Ol = 7] 


Das Symbol ¢ bedeutet den Cauchyschen Hauptwert des Integrals iiber K; das Integral der zweiten 
Zeile stellt den Residuenanteil der Pole bei einem speziellen Umlaufssinn dar, der sich aus asym- 
ptotischen Anforderungen an das Potential ergibt. Der Ausdruck (1) wachst iiber alle Grenzen 
im Falle Q =+ 1/4. 

Die Quellsingularitaten zur Erzeugung der Kérperumstrémung werden in ihrer Ergiebigkeit 
angesetzt in Proportionalitat zur Zeitableitung der értlichen Eintauchung, d.h. zur Strémungs- 
geschwindigkeit normal zur Kérperoberflache. Zur Erfassung der Variation des Deplacements 

fiir die Wellenerzeugung durch Fahrt voraus ist es notwendig, jede solche Quelle zu erganzen durch 
einen Dipol oszillierender Intensitat, dessen Achse in Fahrtrichtung liegt.* 

Ist q ein MaB der ortlichen Tauchamplitude, so wird die Quellergiebigkeit gleich iy q e'”', das 
Dipolmoment gleich — V q e'”‘ komplex angesetzt. Dabei ist y die Frequenz der Schwingung im 
mitgefuhrten Bezugssystem. Das Potential qgq* dieser Quell-Dipol-Kombination erhalten wir 
dann aus @ (1) durch 


ae dp 
ap* =a(iro— eae (2) 


Stellen wir die Zeitableitung im mithewegten System dar durch das Symbol @/dt, die Zeitableitung 
_im ruhenden System durch d/dt, so gilt 


OG) a ; 
ae per Ve (3) 
und damit 
Z 4 
p= 7 (4) 
Zz 3. Energiebetrachtungen. Bewegt sich eine Quelle der Ergiebigkeit iy qe” mit Translations- 


_ geschwindigkeit V im Felde einer mit Frequenz » relativ zur Quelle periodischer Stroémung von 
einem Potential y = ye!”', so ist die aus dem Raum um die Quelle pro Periode abwandernde 


Energie** gegeben durch 


$ .., dw 
Nhu = 4 Re loin ges S, (5a) 


* Siehe FuGBnoten 3 und 4 von Seite 39. 
** Siehe FuBnote 2 von Seite 39. 
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Dem abwandernden Impuls pro Periode entspricht eine Kraft P auf die Quelle in x-Richtung von 


der Grobe 


Tee emie 0p 
pia—1 Reloingetr 2, (5b) 
dem wiederum ein Leistungsanteil 
1 , ay 
ve ae ee Figu yt 5¢ 
Ny 5 Re|gin Ve = (5c) 


entspricht. Der Dampfungsenergie entspricht die Differenz von (5a) und (5c) 


; ., Ow 1 9 — 
Nb =4 Relo in gett F| = 5 Re {oqo “ 


Diese Energie muB je Periode aufgebracht werden, um die Schwingung aufrecht zu erhalten. 
Analog erhalten wir fiir den Dipol der Intensitat q V 


Ne: == > Re l V qei’! i a (6a) 

eee: vig eine (6b) 

Np= B Re lo Voges = ae (6c) 
und damit erhalten wir fiir das Quell-Dipol-System 

Nid == Re loge os (7a) 

Nii = — + Re log ven 2H, (7b) 

Ni’ == Re gg eet a Ht : (7c) 


Ist wy die durch die Singularitaten erzeugte Strémung, so ist in (5) einzusetzen y = q :; in den 
Formeln (7) hingegen y = qp* = q dy/dt. Es ist dann (r Eg) = — 2 f zu setzen. Es leistet dann 
aber nur das Integral in der zweiten Zeile von (1) einen Beitrag zu den Ausdriicken (5), (6) und (7). 
Dieses Integral bestimmt das Wellenfeld in groBer Entfernung. 


Nun ist aber! 


Sb, oni False 5; i (VK; cos @ +9) eit Fo) t ir _ ois Vg K; eXi(r Eo) +ivt Gi as 1:2y5, (8) 


Damit erhalten wir aus (1) und (5) fiir die Quelle in der durch sie erzeugten Strémung vom Poten- 


tial qq 


3 “ | 
Nb= G99 Ke [z= ra Ke em? Sf 6g Ke.) dO, (Qa) a 
; No oF aie sec O 2 > — 2 K 2) ie - 
Wg to So) ag, OE ee ee ee 
ie loa - (6) 
Ngo =" 19 Ko ae K, (Kp? oP SF + KP eo? KS) dO. Coe 
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: Fiir die Quell-Dipol-Kombination ergibt sich hingegen aus (1), (4) und (7) 


q@ eg = sec? O Z . 
Nd = %,Koqav a Ki e* SF 4-6, K2 e—?®&f) dO, (10a) 
0g = é sec O 3 < 
Ni = ex M094 | eo (0, Kp e~? "Ff + 6, K3e &f) dO, (10b) 
d 0 = ‘ 20 - 
Nhs = £8 gl Ky qq | Ro (Ki eS + KEI dO, (10¢) 


—n 


Formel (9a) ergibt mit K,->0o die Dampfungsleistung einer ruhenden oszillierenden Quelle, 
_ der Beitrag von K, verschwindet und es wird K, = K= v?/g. Im Falle y = 0 hingegen ergibt 
(10b) die Widerstandsleistung eines Kérpers, der durch einen diskreten Dipol vom stationaren 
Moment q dargestellt wird, der Anteil von K, verschwindet und es gilt K, = K, sec? @. Im Sektor 
|O| = a/2 leisten die durch K, vertretenen Fahrtwellen einen negativen Beitrag zur Dampfung, 
die durch K, vertretenen Schwingungswellen einen negativen Beitrag zum Widerstand*. 


: Berechnet man die Dampfungsleistung nach Haskind und Havelock durch Druckintegration 

_ des Potentials einer Quellverteilung ohne Dipolzusatzterm, so erhalt man Ausdriicke von der 

_ Form (9c); da die Translationsgeschwindigkeit zwar in der Druckgleichung, nicht aber im Sin- 
gularitatensystem beriicksichtigt ist, liegen diese Ausdriicke in den Potenzen der Wellenzahlen 
zwischen (9a) und (10a). 


4. Fahrt gegen regulares Wellensystem. Nun sei ein Wellenfeld gegeben durch 


EIS ore or (11) 


ai) 


wobei r der Ortsvektor beziiglich eines Koordinatenursprungs auf der freien Oberflache z = 0 
ist; E, sei der Vektor (i cos x, i sina, 1) - « stellt den Winkel der Fortschrittsrichtung der Wellen 
gegen den Schiffskurs, d. h. gegen die x-Richtung dar. 1 ist die Frequenz der Welle in ruhendem 
Bezugssystem, K = 2/g die Wellenzahl und h die Wellenamplitude. 


In einem mit Geschwindigkeit V mitgefiihrten Bezugssystem hat die Welle die Begegnungs- 
frequenz 
: Re|ga 2 Q 12 
y = Re al ) COS &) (12) 
— mit 2, = » V/g. 
| Durch (12) wird der Bewegungsfrequenz ein negatives Vorzeichen zugeordnet fiir den Fall, . 
daB die Anderung des Schwingungsausschlags gegensinnig zu der Anderung der Wellenhéhe im 


ruhenden System erfolgt. 
Z Aus (12) folgt 
. Q cosa = Q, cos a (1 — Q) cos x) . (13) 


Unabhangig von der Fortschrittsrichtung « des Wellenzuges kann also Q cosa nur im Bereich 

0 <Q, cosa < 1 positive Werte annehmen, und zwar den maximalen Wert 1/4 bei Q, cosa = 1/2, 

also wenn das Schiff mit der halben Geschwindigkeit g/y der tberholenden Welle fahrt. Die Auf- 
 lésung von (13) nach Q) cos a ergibt 


; 1 
2 Qo cos 0 = > 


(1 +1 —4.0 cos) ; (14) 


= * Siehe FuBnote 1 von Seite 42. 
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damit wird 


see ae rele 
(1— 2 cosa)? =|3 (14 /1—4Q cosa) (15) 
und somit 
y2 ye 2 2 
Dee a ie er (i 


d.h. mit (1b), (1c) haben wir, da K, und K, durch Vorzeichenwechsel der Wurzel auseinander 
hervorgehen: 

Ko == 1G) (17a) 
bzw. 


K =K,{a). (17b) 


je nachdem, ob 2, cosa gréBer oder kleiner als 1/2 ist. 


Wir wollen nun die Uberlagerung betrachten von einem Wellenfeld (11) mit dem Feld gp baw. 
qp* einer oszillierenden Stérung; insbesondere wollen wir annehmen, daf die Oszillation durch die 
Wirkung des Wellenfeldes (11) stationar aufrecht erhalten wird, daB es sich namlich um das Os- 
zillationswellensystem eines Schiffes im Wellensystem (11) handelt. Dem abwandernden Impuls 
des resultierenden Wellensystems entspricht der Mittelwert der Kraft zur Uberwindung der Wider- 
standserhéhung durch Seegang und Oszillation; der hieraus resultierende Leistungsanteil Ny muB 
aber der gesamten abwandernden Leistung entsprechen, wenn die Energie zur Aufrechterhaltung 
der Schwingung von der Welle selbst hervorgebracht, also nicht von auBen zugefiihrt wird. Es 
muB also Np, die Differenz von Gesamtleistung und Widerstandsleistung, verschwinden, d. h. 
aber, daB der Anteil Np, des Wellenfeldes (11) dem negativen Anteil Np, der Oszillationswellen 
gleich sein mu. Fiir das Potential (11) ist nun aber bei beliebiger Quell- und Dipolverteilung aus 
(5) und (6) ersichtlich, daB zwischen der zugehérigen Dampfungsleistung Np, und der Widerstands- 
leistung Ny, die Beziehung besteht 


1 1 
ip ND. = GV Keosa (18) 
Damit kénnen wir aber die gesamte Widerstandsleistung durch 
V Kos « 
Nyy ges = Ny, — ” ==> ING (19) 


d. h. durch Widerstandsleistung N,,, und Dampfungsleistung Np, des Oszillationsfeldes ausdriicken. 
Kinsetzen von (16) ergibt dann 


2 2 : 
Ni ges = Nir, 2 0080 (p= Np, + (20) 


Im Falle 0 < Q, cosa <1, d. h. wenn das Schiff vor der Welle lauft, aber noch nicht iiberholt 


wird, fallt der zweite Term von (20), welcher der Widerstandserhéhung durch den Seegang ent- _ 


spricht, negativ aus, da in diesem Bereiche Q cosa = 0 ist. Im ibrigen erscheint bei gleichem 


{Q, cos « die Widerstandserhéhung gréer wenn Q, cos x > 1/2 ist, wenn das Schiff also die Wellen — 


einholt, als bei Fahrt gegen die Wellen; im letzteren Fall ist bei gleichem V die zugehérige Fre- 
quenz v, kleiner, d. h. die Wellen langer. 


Es hat den Anschein, als ob in dem Ausdruck (20) die Phasendifferenz zwischen der erregenden 
Welle und der oszillierenden Singularitat nicht eingeht; im Falle eines schwingenden Koérpers ist 


aber diese Differenz gerade bestimmt durch die Dampfung der Schwingung, welche den zweiten 
Summanden von (20) ausmacht. 


Setzen wir den Ausdruck (9) bzw. (10) in (20) ein, so erhalten wir 
Na ov V fe * sec? O 


W ges — 82 Ky K, — K, 


—nx ° 


x [6, K, (K, cos O — K cosa)e 7 ™F 4 6. Ky (K, cos 9 — K cosa) e ***4] dO, (21a) 


a 
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W ges Kok, 


sr 


Ce leae. : sec? O 
Ni = "fe Koa [ 


x [6, K?(K, cos @— K cosa) e °F + Os K? (K, cos 0 — K cos x) e  *S]d0. (21b) 


Pur «x =, Q>0 d =) : oe S. ‘ 
lols a ale oder « , 2 > 0 gilt nun aber: — K cosa = K, (a) bzw. K,() und da fir 


|K, cos O| < K.(z) < K,(2) S |K, cos O| (22) 


so folgt, daB auch in diesem 9-Bereich die Betrage der Summanden von (21) positiv sind, daB also 
der Gesamtwiderstand bei Fahrt senkrecht zur Wellenfront stets positiv ausfallt, solange die 


Welle das Schiff nicht tiberholt. 


5. Widerstandsbeiwert und Dampfungsbeiwert. In den unseren Untersuchungen zu Grunde 
liegenden praktischen Problemen haben die gesuchten GréBen nicht die Dimension von Leistungen. 
Fir den Widerstand interessiert der Mittelwert der Kraft zur Uberwindung des Zusatz-Wider- 
standes, fiir die Betrachtung der Schwingung die Abklingkonstante einer freien Schwingung. 


Wir fuhren dimensionslose Parameter von Fahrt und Frequenz ein durch 
XY == Sun — eee 9 
VKof VfFe 


(23a) 


N=VK f=" Vf - (23b) 


Es gilt dann Q = ( %. Wir haben als Bezugslange fiir die Frequenz die Eintauchtiefe f der Stérung 
gewahlt, weil das Wellenbild, insbesondere das relative Verhalten von Dampfung und Widerstand, 


' entscheidend von K f, d. h. von der Tiefe der Singularitaten, abhangt. Auch bei Schiffen erscheint 


es deshalb angebracht, die Frequenz nicht durch die Breite**, sondern durch den Tiefgang dimen- 
sionslos darzustellen. 


Einen dimensionslosen Dampfungsbeiwert definieren wir dann durch 


_ 8x fi? 1 
Des agate ye ND 9 (24a) 


einen Widerstandsbeiwert W, der erzwungenen Schwingung durch 


8a fs? 1 
M1 = ogg? BM (24b) 


- und den Gesamtwiderstandsbeiwert W gem4B (20) durch 


2 
I= J1—42 cos « 


W=W, =A ) MN? D, cosa. (24) 


~ Die Formeln (5a), (6a) zur Ermittlung abwandernder Energie gelten in verallgemeinerter Form 


_ dann, wenn die Kérper nicht voll eingetaucht sind, sondern die freie Oberflache durchstoBen 


fiir alle bewegten Kérper, die durch Singularitatensysteme dargestellt werden, insbesondere auch 
*** 


Es braucht dabei nicht die Annahme eines kleinen Breiten—Langen-Verhdltnisses gemacht werden, 


_ solange sich die entstehenden Wellenfelder nur linear superponieren lassen. Der Ausdruck fiir die 


gesamte Widerstandserhéhung W im Seegang, der aus dem Ansatz (5a) baw. (5b) berechnet ist, 


ee gilt damit unter denselben Bedingungen auch fiir nicht voll eingetauchte Kérper. 


* Siehe FuBnote 1 von Seite 42. 
** Siehe FuBnoten 1 von Seite 39 und 3 von Seite 40. 
*** Siehe FuBnote 2 von Seite 39. 
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Fiihren wir noch dimensionslose Wellenzahlen ein durch 


@ = Kh (25a) 
yes (25) 
Urea 6) fe (25c) 
‘so erhalten wir aus (9) und (10)-far die Quelle 
D,= | gma HT $82) do , (26a) 
t J 2 
= 5 = (5, RZ e- 2% + dy KZe-2%) dO, (26b) 


Y= ee | ae [d, R, (RK, cos O — & cos x) e— 2% + 6, Ky (Ky cosO — Reosa)e—?%] dO (26c) 
1 2 


HS 


und fiir die Quell-Dipol-Kombination 


aioe J sec? oO (5; RZ e— 2% + 5, F2 e— 2&2) dO (27a) | 
q MB? ee Bi 
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Abb. la und b. Dimpfungsbeiwert D, in Abhangigkeit von Frequenz und Fahrt, — a) fiir Quelle; b) fiir Quell-Dipol-Kombination. tt | 
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he 
05 


Abb. 2a und b. Widerstandsbeiwert W, der erzwungenen Schwingung in Abhangigkeit von Frequenz und Fahrt. 
a) fiir Quelle; b) fiir Quell-Dipol-Kombination. 


BS 
We, = 5 [ gm, (51 Met + dy Me?) AO, (tb) 


eth 2 


W = a5 | = a [ds HFG, cos —Reosa)e~*% + 5 (Ms 0086 — 8 cosa) e~ #8]. 40 , (27c) - 


at 


Die Ausdriicke, welche Haskind fir D, und W, durch Druckintegration des Potentials der Quelle 
erhalt, liegen in der Potenz der Wellenzahlen zwischen (26ab) und (27ab) in Analogie zu (9c). 


Die Ausdriicke (26) und (27) wurden fiir eine Reihe von Wertepaaren Jt, % berechnet, wobei 
Q die Werte 0; 0,5; 1; 1,5; 2,5; 5 und K, f die Werte 0,1; 0,2; 0,3; 0,6; 0,9 undoo durchlief. Die 
Integrale wurden nach der Simpsonregel mit Schrittweite 7/18 berechnet, wobei die Umgebung 
Jer Stelle 2 cos @ = 1/4 besonders abgeschatzt wurde. Der Ausdruck fiir W wurde fir den Fall 
berechnet, daB die Fahrt gegen die Wellen erfolgt (x =, © = 8,(z)). 


Die errechneten Werte fir D, nach (26a) und (27a) sind in den Abb. la und 1b aufgetragen 
iber N. Obwohl fiir Q ¥% — 1/4 die Ausdriicke fiir D und W iber alle Grenzen wachsen, lassen 
ich die Punkte glatt verbinden mit den Werten fiir J{ = 0; der Bereich der Spitzen, in denen das 
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_Anwachsen erfolgt, erscheint sehr schmal, wie auch Havelock genauer festgestellt hat. Eine Auf- 
tragung von N> itiber & mit zusatzlicher Berechnung der Werte fiir N & — 0,245 ergab noch eine 


monotone Tendenz der Kurve; fiir Ny hingegen lieBen sich die Punkte nicht durch glatte Kurven 
verbinden! 


Abb. lb zeigt, daB der zusatzliche Dipolterm in (27a) fast im ganzen Bereich (bis auf Q = 0,5; 
K,f = 90,9 und 0,6) eine Verminderung der Dampfung mit zunehmender Fabrt gegeniiber den 
Werten fiir die Quelle nach (26a) bringt; schon aus diesem orientierenden Vergleich kénnen wir 
den Schlu8B ziehen, daB der Ansatz von Haskind das Abklingen der Dampfung mit zunehmender 
Fahrt zu schwach erfaBt. 


Dem starkeren Abklingen der Dampfung beim Quell-Dipol entspricht eine VergréB8erung des 
_ Widerstandsbeiwertes W, fiir die erzwungene Schwingung; (der Anteil der Fahrtwellen gewinnt 


durch den Faktor §,/Q® starker als die Schwingungswellen durch den Faktor Ks /N?). 


Dagegen klingt der zweite Widerstandsanteil, W— W,, proportional zur Dampfung starker 
ab als bei der Einzelquelle. 


6. Darstellung schwingender Kérper durch Singularititenverteilungen. Wir dehnen unsere 
Untersuchung aus auf kontinuierliche Verteilungen, welche sich langs einer horizontalen Geraden 
—Lj2sx*xs 1/2, y=0, zs =—f erstrecken und die Umstrémung eines Schiffskérpers der 
Lange L approximieren sollen. Fiir eine Quellverteilung der Ergiebigkeit iy q e’”‘ o(x)/L erhalten 
wir dann nach Aufstellung des Geschwindigkeitspotentials und Berechnung der Leistungsanteile 
-unter Beriicksichtigung des Dipolzusatzes 


P24 


1 260 3 a 
D, = gus | C=, (6, 93 HH, + 6, 93 H, H,) a0 (28a) 
| und 
1 O =e =e 
W= gs | GES 98 HH + 85.9) HH) 40 (28) 
mit 
1 
Bess eels et ane de (j= 1,2), (29) 


—l1 


wobei € =22/L, 1 = L/2 f ist. 


Fiir Tauchschwingungen eines quaderférmigen Korpers mit kleiner Ausdehnung in y-Richtung 
setzt man an o(x) = 07(x) = 1; es wird dann 


H, = e— 8 Ny (8; | cos O) (Gia 1-2) (30a) © 
| mit 
% 1 ce 
: N,(a) = | cos (Rf) dé =— “2. (30b) 


0 


Fiir Stampfschwingungen desselben Kérpers, welche den Tauchschwingungen um € — 7/2 in der 
| Phase vorauslaufen mégen, wird analog gesetzt 


o(x) =0,(x) =e °—™ ; (31a) 
es wird dann 


Hj =e %i +** M, (®; 1 cos @) (j = 1,2) (31b) 
J 1 (Sj 


PEER BOR ai r 
| 
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mit 


M,(A) = f €sin (A) dé. (31¢) | : 


Fir diese beiden Falle wurden die Gré8en D,, W, und W nach (28) und (24c) berechnet. Bei i 
Uberlagerung von Tauch- und Stampfschwingungen sind die entsprechenden Verteilungen a(x) | 
zur Bestimmung von (29) — eventuell mit verschiedenen Gewichtsfaktoren— zu superponieren ; 
in den Ausdricken fiir D,, W, und W treten dann noch gemischte Glieder auf, welche den Faktor 
cos € enthalten, also verschwinden, wenn die Tauchschwingung mit der Stampfschwingung syn- | 
chron ihre maximalen Werte erreicht. (Haskind betrachtete in seiner Analysis nur diesen Fall } 


D 
7 


Abb. 4a —c. Dampfungsbeiwert des schmalen Quaders in Abhiangig- 
keit von Frequenz und Fahrt. — a) Tauchschwingung; b) Stampf- 
Abb, 4c, schwingung; c) Kopplung. 


und vernachlassigte deshalb in seinen Bewegungsgleichungen den sich aus D, ergebenden Beitrag 
zur Dampfungskopplung.) Die Kopplungsterme lauten 


ja fede ZO | 
sas ne | ®, = @, (1 Bi No My e~*® + 52 8; No My e~**) dO, (32am : 


1 a 
Wi =e i a (0, &t No My e~? % + 6, 83 Ny M, e—?%) dO. (32b) — 


—T 


Die Ausdriicke (32ab) verschwinden (und zwar allgemein fiir in Langsrichtung symmetrische | 
Schiffe) identisch fir 8 —=0. Der Audsruck D, cose gibt nur den Teil der Geschwindigkeits- — 
kopplung der Differentialgleichungen fiir die Schwingung wieder, welcher zu Energiedissipation 4 
fiihrt, der also in beiden Gleichungen dasselbe Vorzeichen hat. z 


Fir die Singularitatenverteilung (31a) hat Havelock fiir 1 = 5 den Ausdruck fir D, berechaea 
unter Auslassung der Dipolkorrektur. Seine Untersuchungen gelten vornehmlich dem Anwachsen 4 
des Dampfungskoeffizienten im Bereich Q ~ 1/4. Havelock wies darauf hin, daB D, (und damit — 


> 


vin 
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auch W) nur dann fiir 2 = 1/4 unendlich groB wird, wenn nicht M, (K;(0) 2) verschwindet im Falle 
der Stampfschwingung ; fir die Tauchschwingung bleibt D endlich, eee Np (K;(0) 1) verschwindet 

Die folgenden Untersuchungen behandeln den Fall 1 = 10 unter Berick eae der Di Be 
korrektur. Die Rechnungen wurden beschrankt auf Wertepaare MN, %, welche bei ‘Fahrt aide 


Abb. 5a. 
{ a 2 
°B= 316 
9g 0,24 
a 183 
a 409 
e 105 
0 ees eek 
== Z 3 
-01 JE 
Abb. 5b. 
02 ad Abb. 5a — ec. Widerstandsbeiwert der erzwungenen Schwingung des schmalen Qua- 
ders in Abhangigkeit von Frequenz und Fahrt. — a) Tauchschwingung; b) Stampf- 


Abb. 5c. schwingung; c) Kopplung. 


Die errechneten Dampfungsbeiwerte D,, nach Formel (28a) bzw. (32a) unter Benutzung von 
(30a) und (31b), sind in Abb. 4a, b, c dargestellt; die Kurven fiir %§ = 0 sind in Abb.7 noch einmal 
wiederholt in Gegeniiberstellung ‘zu der entsprechenden fiir die Quelle. Abb. 7 laBt erkennen, 


_ daB®B die gewahlte Langenausdehnung der Quellverteilung im Vergleich zu den Wellenlangen, welche 
_ zur Dampfung wesentlich beitragen, so klein ist, daB die Kurven fiir Quelle und Tauchschwingung 


weitgehend iibereinstimmen; deshalb erscheint auch der Dampfungsbeiwert der Stampfschwingung 


in diesem Bild von wesentlich kleinerer GréBenordnung. Die ausstrahlenden Wellen sind hier 


durch Interferenzen abgeschwacht, welche fiir die Tauchschwingung ‘erst durch den FahrteinfluB 


: _ wirksam werden. Der Kopplungsanteil, welcher fiir 3 = 0 identisch verschwindet, bleibt auch in 


Fahrt unbedeutend gegen den Anteil der Tauchschwingung. 


Die entsprechenden Widerstandsbeiwerte W, der erzwungenen Schwingung, berechnet nach 


- Formel (28b) bzw. (32b), sind in Abb. 5a, b,c dargestellt. Hier erscheint der Koeffizient der 


b> 


_ Kopplung dem der Tauchschwingung in der GréBe vergleichbar. 


Die Widerstandsbeiwerte W der freien Schwingung in der Welle finden wir fiir den Quader 


auf Abb. 6a, b, c. Wiederum hat der Kopplungsterm die GréBenordnung des Tauchschwingungs- 


Seer) ON Le eee ere. 
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Abb. 6c. Abb. 7. Dampfungsbeiwert D, in ruhendem Wasser (% = 0) in Ab- 
hangigkeit vom Frequenzparameter Jt. — + Einzelquelle (Unterer 
Teil der Kurve enthalten im Abb. la und 1b); © Tauchschwingung 
Abb, 6a —c. Widerstandsbeiwert der freien Schwingung des schma- des Quaders (Unterer Teil der Kurve enthalten in Abb. 4a); 
len Quaders in der Welle in Abhangigkeit von Frequenz und Fahrt. * Stampfschwingung des Quaders (Kurve in gréBerem MaBstab 


a) Tauchschwingung; b) Stampfschwingung; c) Kopplung. enthalten in Abb, 4b). 
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Abb. 8. Ubersicht des untersuchten Frequenz- und Fahrthereichs. 


(A = 25f) Bereich, in dem Haskind fiir Michell-Schiffe D, und W, berechnete unter Vernachlassigung der Dipolkorrektur*. 


(A = 32f) Bereich, in dem Maruo die gesamte Widerstandserhéhung bei Fahrt gegen Wellenzug vorgegebener Héhe und da- 
durch angeregte Tauch- und Stampfschwingungen fiir eine Schiffsform berechnete, die er durch oszillierende Quellen in den 
Schwerpunkten der Spanten darstellte (** und Int. Shipbuilding Progress. 4 (1957) S. 477). 


(2 = 0,66) Bereich, in dem Hanaoka fir ein Michell-Schiff in erzwungenen Stampfschwingungen D, und W, berechnete. 
(Vorgetragen auf Tankleiterkonferenz Oslo 1954.) 


Bereiche, in denen Hanaoka experimentelle Messungen des Gesamtwiderstandes im Seegang vornahm. Die zugehérigen Be- 
rechnungen basieren auf der Annahme, daB der Beitrag des Seegangs zum Widerstand durch die Annahme einer wellen- 
férmigen Verformung des Schiffskérpers erfaBt werden kann. Es ist hier f = L/25 angesetzt***. 


Bereiche, in denen Havelock fiir die Stampfschwingung (Singularitatenansatz nach unserer Formel (3la)) unter Vernach- 
lassigung der Dipolkorrektur D, berechnete****. 


Bereiche, in denen Golovato ***** fiir ein erzwungen stampfendes Modell D, experimentell ermittelte. 
Punkte, fiir die in der vorliegenden Arbeit W,, D, und W berechnet wurden. 


Errechnete Grenzen des Bereichs negativer Widerstandserhohung W, in erzwungener Schwingung. I fiir Quelle; II fiir Quell- 
Dipol-Kombination; III fiir Stampfschwingung. 


i Symbole: 9 =v /f/g = Frequenzparameter 


&@ = V//'gf = Fahrtparameter 
2 =NB = Brardscher Parameter 


* S. FuBnote S. 40 

** S. FuBnote S. 40 
*** S. FuBnote S. 39 
*#** S. FuBnote S. 39 
e#4##* S, FuBnote S. 40 


f = Tauchtiefe der Singularitaten bzw. Modelltiefgang. 
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terms, diesmal mit definit positivem Vorzeichen. Die Kurven fir ¥% = 0, welche in Abb. 6a, b 
3a, b nicht voll ausgezeichnet sind, ergeben sich aus den Kurven in Abb. 7 nach Multiplikation der 


Ordinatenwerte mit Jt? gemaB Formel (24c). 


In Abb. 8 finden wir zum AbschluB einen Uberblick iiber die Bereiche von %t{ und %§, in denen 
bisher rechnerische und experimentelle Untersuchungen tiber Dampfung und Widerstandserhéhung 


bekannt sind. 


Die umfangreichen numerischen Berechnungen obiger Ergebnisse wurden von den Herren 
cand. Jeremie, Béhm und Kleuters am Institut fiir Schiffbau, Hamburg, durchgefiihrt, erméglicht 
durch eine Zuwendung der Deutschen Forschungsgemeinschaft. 


Der Verfasser méchte hier seinen Dank aussprechen fir diese sachliche Hilfe, ebenso jedoch 
an Herrn Prof. Dr. Weinblum fiir Anregung und Férderung obiger Untersuchungen. 


(Eingegangen am 23, April 1959.) 


Anschrift des Verfassers: Dr. rer. nat. Dipl.-Math. Klaus Eggers, Hamburg-Volksdorf, Horstlooge 16. 
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The plastic yielding of double notched bars due to pure bending 
By George Lianis 


1. Introduction. In a previous paper the author! has treated both theoretically and experi- 
mentally the problem of plastic bending of single notched bars. In the case of double notches, 
when pure bending is applied, the stress pattern and the plastic nucleii developed under the notch 
roots are exactly similar, with opposite signs, since, due to the symmetry, the upper half of the 
bar suffers as much in tension as the lower half suffers in compression. The two similar nucleii 
developed during the elastic-plastic stage of bending spread into the elastic material simultaneously 
with increasing bending until they join together. It is this value of the bending moment which 
is defined as the Yield or Collapse Moment. 

For the yield moment a slip-line field is constructed for any shape of symmetrical notch. It is 
proved that this field does not involve anywhere negative rate of plastic work. A velocity field 
associated with this slip-line field can be constructed. The latter is consistent with the boundary 
kinematic conditions and the rigid-body motion of the assumed rigid parts. The moment calcu- 
lated from the associated stress field is strictly an upper bound.? A suitable stress field satisfiying 
static conditions is also constructed by means of which a lower bound is calculated. The special 
cases of Rectangular, Trapezoidal, and V-notches with circular fillets’are examined completely 
and lower bounds are found for V- and Circular notches for which Green* has found an upper bound 
solution. 

In a series of experiments for the above shapes and various sizes it has been found that the 
measured constraint factor agrees satisfactorily with the upper bound. 


2. Limit analysis. We are concerned here with initial motion problems only. It is therefore 
permissible to satisfy the boundary conditions on the undeformed boundary. A complete solution 
for a plane strain problem must satisfy the following conditions: 


a) The equilibrium equations 


00% Oxy aoy Oty 
at Saas and rope oe == 0. (1) 


b) The yield criterion throughout the plastic zone* 
(6, —o,)? +472, = 41 (2) 
where k = Y//3 , and Y = the tensile yield stress. 


c) The stress system in the plastic region must not cause violation of the yield criterion, equ. (2), 
in the rigid parts. 


d) The flow rule must be satisfied by the strain system 
P 


Se = ey = Vey (3) 
Ox—Oy Tay 
and the equation of incompressibility 
P+ &P — 0 (4) 


in the plastic region. 
e) The velocity field must determine in the rigid region rigid body motion. 
f) The rate of plastic work must be positive everywhere in the plastic region. 
g) Stress and velocity boundary conditions must be satisfied. 


Any solution which does not satisfy but some of the above conditions is called incomplete 
solution. Incomplete solutions are sometimes useful to obtain bounds to the yield load. 


1 G, Lianis and H. Ford, J. Mech. Phys. Solids 7(1958) p. 1. 

aD): C. ee H. J. Greenberg, E. H. Lee and W. Prager, Proc. Ist U.S. Nat. Congr. Appl. Mech. 1952. 
3 A. P. Green, Quart. J. Mech. Appl. Math. 6(1953) p. 223. 

4 Von R. Mises, Nachr. Ges. Wiss. Gottingen, math.-phys. Klasse (1913) p. 582. 
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We shall call a stress field statically admissible if it satisfies conditions a, b, c and g, and || 


a velocity field kinematically admissible if it satisfies conditions d, e, f and g- ; 
The following theorems have been formulated for the case of all surface tractions increasing 
proportionally?: me: 
Theorem I: The yield load is the largest between the loads for which a statically admissible 
field can be found. 


Theorem 2: The yield load is the smallest between the loads for which a kinematically ad- — 


missible field can be found. 


A statically admissible stress field is described here by means of which a lower bound for the 
yield moment is found. The slip-line field theory on the other hand provides kinematically ad- 
missible velocity fields. In a material stressed plastically under conditions of plane strain the 
velocity and stress characteristics coincide at each point of the deformed region with the directions 
of the maximum shear stress (or shear strain). The two orthogonal families of characteristics in 
this case are called slip-lines. Let the two families be labelled by the parameters a and f 
respectively* (Fig. 1). 


Fig. 1. Fig. 2. 


In a small element bounded by two pairs of «- and f-lines, p = = (0, + 6.) is the normal 


stress on each face known as hydrostatic pressure, and k is the shear stress where k = Y/V3, 


Y = the yield stress in simple tension. The following equations known as Henckys equations? — 


relate p with the angle @ of the slip-lines with a fixed system of axes of reference: 
p+2ko = const. along an a-line , 
p—2kqg = const. along a f-line. 


These equations are equivalent to the equations of equilibrium. 


The velocity components u and v in the a- and f-directions respectively are related with the — 


angle » by Geiringers equations! 
du —-v dp = 0 along an a-line , 
dv + udp = 0 along a f-line . 


of the point P (Fig. 1) referred to axes passing though some fixed point 0 and parallel to the 


1 See footnote 2 page 55. 
* R. Hill, The mathematical Theory of Plasticity, p. 134. Oxford 1950. 


* If the a- and £-directions are regarded as right-handed axes of reference the line of action of the alge- 
braically greatest principal atress (0, or 62) lies in the first and third quadrants. 
° H. Hencky, Z. angew. Math. Mech, 3(1923) p. 241. 


“ H. Geiringer, Proc. 3rd Int. Congr. Appl. Mech., Stockholm 1930, vol 2, p. 185. 


(5) 


(6) 


A point in the slip-line field can be defined by the quantities x and y. They are the co-ordinates 
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respective slip-line directions at P (Fig. 1). The variation of % and y along the sliplines is given 
by the equations! 


(7) 


A solution based on a slip-line analysis satisfies the conditions a, b, d, e and g. Ifthe rate of plastic 
work is everywhere positive the slipline field provides a kinematically admissible velocity field. 
The load calculated from the associated stress field is therefore an upper bound to the yield load. 


3. General solution. a) Upper bound. 1.) Slip-line-field; Solution I. The construction of the 
slip-line field is based on the same method followed in the general case for single notches?. The 
notches are similar and symmetrical with respect to the bar axis y’ — y (Fig. 2). Pure bending is 
applied at the ends of the bar tend to open the upper notch. The stresses on the upper half will 
therefore be tensile and these of the lower half will be compressive. Apart from the sign both 
configurations must be similar due to the symmetry. The axis of the bar is therefore axis of Sym- 


dy + x dy = 0 along an a-line, 
dx — y dy = 0 along a f-line . 


_ metry of the slip-line field. 


The free surface of the notch defines as previously was sought the field in the domain 
ZPLP’2Z'L' F’NFLZ (Fig. 2). The field is extended until it reaches the bar axis y’ — y from 
ae sides. The axis of the bar defines therefore the position of the end of the field Z at the notch 
surface. 

The equilibrium conditions are identically satisfied in such symmetrical construction. The 
point NV is a point of stress discontinuity. The validity of this solution depends on whether the yield 
criterion is violated at the corner of the rigid part of the bar near N. If @ is the angle of a generic 
point on the notch surface with respect to y-axis, the slip-line angle on the notch will be 


It 
OF 


Applying Henckys theorem, equ. (5), along ZN 


—k—-2k (02+ 7)=py—2kZ, 
4 
Py =—k(1+207). 
In the compressive part on the other hand it is 
Pyr=+k(1+20,). 


The yield criterion is not violated at N if® 
|Ap| =2k(1+20,)Skax 


or 


0, <7>* = 16°’. (8) 


The limiting values, of a, a);,,3;. corresponds to 97 = 16° 20’. Plotting the normal stress o, applied 
across 0A N and taking the moments with respect to N the yield moment My is found 
1/2 


a Oye | 2) 2 
My =2{ 3 Glan 
0 


The Constraint Factor is defined in the following as the ratio of the yield moment of the notched 
bar to the yield moment of a uniform bar having the same thickness. The later according to the 


simple beam theory is equal to ka*/2. Therefore: 


1/2 


bape a4 Pee(G): (9) 


~ Kk at k a 
hy 


where x the distance measured from the point N (Fig. 2). 
2. Slip-line field, Solution II. If extending the field to reach the middle point N in Fig. 2 a point 


ma — 2 


4 


Z is found on the notch surface such that 0, > in solution I, the yield criterion is violated 


1 See footnote 2 page 56. 
2 See footnote 1 page 55. 
3 R. Hill, J. Mech. Phys. Solids 2(1954) p. 278. 
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in. some parts of the assumed rigid material. Another solution satisfiying the stress and velocity 
boundary conditions and holding for values of the minimum bar thickness greater than @j,,,j; 18 
shown in Fig. 3. This solution yields lower values for the constraint factor than solution [. 

The field is extended as before until the point Z having 9, = (7 — 2)/4 is reached and extended 
further to the left and to the right as before (regions Z TEN, Z’ T’ E’ Nin Fig. 2). The two domains 
joint together by the circular ares EE, and E’ EF’, forming a central hinge. These arcs have radius R 
and centres K, and K, lying on the axis of symmetry of the bar y’ — y. The stresses on points P, P,, 
lying on the boundary of the hinge symmetrically to y — y’ have equal and opposite projections 
in the y-direction and therefore the equilibrium condition is identically satisfield. 


(Il) 


Fig. 3. 


The second equilibrium condition (in the x-direction) along ON EE, N, O, and the fact that K, K 
lie on the y’ — y axis define the two unknowns of the problem, namely the hinge radius R and the 
position of the end points T, T’, T,, Tj. Analytical expresions to these two requirements will be 
given in the special cases examined in the next paragraphs. 


Once the field has been fixed the constraint factor can be evaluated by taking the moments 
N E 
of the stresses distributed along ON (¥) Pe INGE: (%) ,and ES (1) 
to) N E 
N E s 
4 (3 ee Mm) 
0 N E 
k a ; 


The hydrostatic pressure at the points S and S, is zero as one can see b lines Fi 
tion (5) along ZN and NS 7 ee by applying Henckys equa- 


7) TM 
k= 2k/ ~ + Z)=Pw—2k a, Py =—k= 


i 


(10) 


k= +h > =ps+2hx0, 
and therefore 
Ps =0. 


Fe eee eet eee 


Bi ES 
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! 3) Velocity field (Fig. 4). Let us examine the velocity field (Hodograph) for solution II. The 
rigid parts rotate around K, and K, respectively with an angular velocity w. For simplicity we 
take ® = 1 without loss of generality. If now a velocity field associated with the slip-line field 
is found such that it satisfies the boundary conditions and does not involve negative plastic work 
anywhere in the deforned region (condition f) the load given by the slip-line feline an upper bound 
for the yield load. Such a field is shown in Fig. 4. The Hodograph for points of the regions I and II 
along E'S and E’S’ consists of circular arcs of radius R (since w = 1) with the ae as a centre 
and angular span 4. Along ELT: If K, is taken as the origin then due to the continuity of the 
normal velocity along ELT and the fact that w = 1 it follows from Figs. 3 und 4b: u = y and 


therefore from equs. (6) and (7) v =—%. Therefore the image of EL T in the Hodograph plane 
is a similar curve rotated through 90°, x 


Fig. 4. 


Along NE: From the continuity of the normal velocity it follows v = 0. Applying Geiringers 
| equ. 6 we find along NE, and NE’ u = const. = ug = R. The jump of the tangential velocity 
_ along these lines must remain constantly equal to R!, their images in the Hodograph plane being 
circular arcs of angular span A — 2/4. Once the normal velocities along NE and EL are known 
the Hodograph in NELZ can readily be constructed!. Having the distribution of the normal 
_ components of the velocity along ZL and ZT from the above construction the Hodograph in the 
_ region ZLT can be drawn. On either side of NZ and NZ’ thereis a tangential descontinuity of 
amount R since u = 0 at My and u = R at Nyy. (Au = R and this jump must remain constant 
_ along NZ). (NZ)y; and (NZ’)y being known in the Hodograph, the velocity field in (NZZ’)y1 
_ ean be constructed and extended to meet the boundary ZOZ’. 

- The Hodograph of the lower part of the field is exactly similar to that of the upper part. We 
' shall now give a proof of the following general statement: 

_ For a concave upwards notch having monotonically increasing tangent di- 
' rection (with possible discontinuities) the stress and the velocity system of Figs. 
| 3 and 4 cause nowhere in the plastic region negative rate of plastic work. 

| If s, and sp are the legth of the arc element along «- and f-lines respectively the corresponding 
| curvatures 1/R and 1/S are given by 

we, es (11) 


Es oa ae ee S apse 


1 See footnote 1 page 55. 
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and the rate of plastic work per unit of volume is 


; du ov v u : 
w= Sey Urbeg tae ane () 


Green! has shown that if R’ and S’ are the radii of curvature of the images of w- and f-lines in the 


Hodograph plane 

: S’ R’ 

= (5 —R)- (13) 
If by examining the form of the slip-line field and the velocity field S’/S and —R’/R are found to 
be of the same sign, the sign of w may be found without any laborious calculation of the curvature 
using equ. (13). In Fig. 3 the angle ® increases everywhere in the positive f-direction. 
D 


Fig. 5. Fig. 6. 


Because in the x — y system shown in Fig.5 ®, — ©, = ®, — ®; and due to the fact that 
the boundary 4—3 is principal direction trajectory: ®, = 1/2 (O, + ®,). Hense: ®, — OD, = 1/2 


(,— ®,). But on the notch surface ® =7/4 +O and therefore: ®,— O, = — (0, —O,), 
Since @ is increasing it follows that ©, > ®,. Consequently: 


1 a®D 
SS ae eee 
Saas. () 
For the same reason @ increases along a-lines. Therefore: 
1 a®@ 
R~ 0, =": (b) 


Since in both lines @ increases they are both concave upwards. We have seen that the boundaries — 
ELT (f-line) and E’ L’ T’ (a-line) have similar images in the Hodograph plane rotated through 90° 
in the direction of rotation of the relevant rigid parts of the bar. Their images are therefore concave 
to the right and to the left respectively, and so are all the other lines of the two domains. Therefore — 
the curvatures of the Hodograph are : 


a 
wD 
us 

N 


R’ — asi, < 0 ? (c) 
iL aD’ 
Shen hogs =0, (dg 
From (a), (b), (c) and (d) it follows 5 
So y 
i >0 and — - =a 0 


1 See footnote 3 page 55. 
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Therefore everywhere in the plastic region y => 0 and w=0. The sign of the work of the 
instantaneous shear done as an element crosses a velocity discontinuity must be found separately 
using equ. (12). Since as we cross NZ moving in the positive «-direction the jump of v is positive: 

Av=+R, “ = + 0o 
So, 
and therefore @ = + and ¢0 it is along NE’, NE;, and NEj. 

The above proof holds for solution I as well, where the points N, EF, S, Ey, Ej, Sj and E’ coincide 

with the neutral point N. The field remains unaltered in the regions NZZ’ and NZ,Z;, N in the 
_ Hodograph coincides with the pole and the discontinuity R disappears. 

b) Lower bound. 1) Siress field. Consider a trapezoidal stress field (Fig. 6). Equilibrium 

along the discontinuity AC leads to 


Pi = k [sin? x, tan (x, + 8) + 2 cos? a,], (14) 
sin 2 
= J = : 15 
Pe een (13) 
Equilibrium along CD leads to 
— > — 04> P3 = 2k sino, (sin a, + cos a). (16) 


| Ifwe superimpose the field in region III with p, = 0 (tension) on the field in region IV with p, < 0 
_ (compression) the resultant field is one of pure shear (Fig. 6): 


sin 2 a, sin 2 6 


Sa aad ae err EET i (17) 
The yield criterion is not violated in this region if 
sin 2 a, sin 2 6 at (18) 


sin 2 (@,-+ Bp) 


~ (2ke-Dyrp, sin’) 


IN 


| Ort cos‘) 


“Let us now combine the field of Fig. 6 with a similar field having p, = compressive as shown in 
| Fig. 7a. In the region IV the yield criterion is not violated if equ. (16) holds. The stress com- 
: ponents in region V are indicated in Fig. 7a. The yield condition is not violated in this region if 


2 
| ie k— pz cos 2B)? + ©? cos* 2B aaa 


Aa) B S P, sin 2 a, (19) 


io 
Gy 
b 
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and this becomes maximum for 
B= 15°, «0 = 31.555 slime nae = 1:20). 
Because the upper bound is for this case 1.38 the mean value of the constraint factor will be 
equal to 1.32 with +4.5% error. The present analysis will be employed to obtain bounds for 
various shapes of notches. Inequalities (15), (18), (19) together with equations (14) and (17) are 


simultaneously satisfied if 


eer AG m% 
ee 9 Gi ee 31 9 ad ere Be 
The constraint factor has maximum value when 
M4 . 
Palos, pen wae ee 20) 
From Fig. 7a and equ. (14) and (20) it follows : 
Pi =2k (1 +sinf), (21) 
2 

Lower = ae eer sin /) > (22) : 

B < 15° has a value to maximise L,,,,,,, and it is calculated from the equation 

“(@) 
7 0088 + Gz" (I + sin f)= 0. (23) 


4, Rectangularly notched bars. a) Upperbound, Solution I. 1) Slip-line field. Since the con- 
figuration is symmetrical with respect to the central axis it is sufficient to examine one half of the 
field (Fig. 8). The field is extended around the singularities Z, Z’ by circular fans of radius b /2. 
The field in the region ABN B’ is constructed easily. The angle (f) is defined by the condition 
that the field is completed on the axis of the bar y — y’. 

From equ. (1) it follows 


am —2 ° , 


By drawing the field for various values of f between 0 and (a — 2)/8, c/a can be found. This solution 
is valid for 
f<@, and therefore = = 0 
2) Constraint factor. For each value of 2b/a and therefore of f the stress distribution along AN 


is calculated by means of Henckys theorem. Taking moments of the stresses along A N with respect 
to N the values of the constraint factor have been calculated and plotted in Fig. 10. 


Fig. 9. 


b) Upper bound, Solution II. 1) Slip-line field. In Fig.9 one half of the slip-line field | 


representing this solution is sketched. The angle ©, remains constant for the various ratios 2b/a 


given by equ. (24). This is a consequence of Henckys theorem and the fact that the hydrostatic — 


pressure at S is zero. 


: 
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The equilibrium condition in the y-direction is identically satisfied due to the symmetry. The 
equilibrium condition along OANDS in the x-direction leads to 


D 
Ss 
= leat) 
EO ize I N . 
ae (F-2 2 @,)cos ®,—(% +142 6,)sin 6] = (25) 
Since K, and K, lie on the axis Y’ — Y it follows from Fig. 7 
R IG 1 b x b 
4 8 (¢— 1] = 7-135 — 5 (By as (26) 


where x(®,) is evaluated by drawing the slip-line field. From equations (25) and (26) the two 
unknowns of the field, R/a and ®, in terms of 2b/a can be found. Since 
D 


ges 


for any value of ®, the L. H. S. of equ. (25) must be positive. Therefore: 


5 —1+2o, 

tan D, < ee (27) 
wt 
+ +1429, 


3 0< 6, <21°50'. 
For ®, = 21° 50’, 2 b = 0 and the notch becomes a slit. 


2) Constraint factor. By taking the moments of the stresses distributed along 04 NDS with 


respect to K, and by means of equ. (3) the constraint factor has been calculated and plotted in 
Fig. 10. 


13 


Ei 


o Copper 
4 Aluminium 


LZ 


IE 


Fig. 11. 


c) Lower bound. The stress field of Fig. 5a fit to a rectangular notch is shown in Fig. 11 


pea 
So ; tanp. 


Rs The lower bound for the constraint factor takes the form 


re fe tan B) (1 + sin f) (27) 
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f is found from equ. (13) which takes the form 
cos 8 + cos f sinh + tan B = oi (28) 


and 
Be 


Both bounds with a mean value curve are plotted in Fig. 10. 


5. Trapezoidally notched bars. a) Upper bound, Solution I. 1) Slip-line field. N is the 
neutral point where a stress discontinuity occurs. Again: 


na—2 
as 
By means of Henckys theorem along A BC and CM PL we find 
= GING 
eee | (29) 
Therefore solution I holds for 
Bt t<cyss sor, 1s lee (30) 


The condition that d => 0 further restricts y for each value of 2b/a. Since c/a is a function of f 
and therefore of y, d/a can be found from Fig. 12 


die 1 b b 
(2 =5-7— Pe. (31) 
For dja= 0 
2b 1 
= : (32) 
a “(f) 
14" 


2) Constraint factor. Taking the moments of the stresses along OA DN with respect to N and © 


integrating: (Fig. 12) 
Sie 
2M Z 
ee A ee | ee (33) 


dja V2 


Fig. 12. Fig. 13. 


b) Upper bound, Solution II. 1) Slip-line field. When y S 13° 30’ a central hinge is de- 


veloped. The upper and lower region are similar having their centres on the bar axis. Applying 


He ; ce te 5 ; 
ee nee the angle @, is found: ®, = (7 — 2)/8. This value gives zero hydrostatic pres- 
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Applying Henckys theorem along A BCE: 


eh ate 
Drie ea ellie (34) 
: 


®, = 


Due to the symmetry the equilibrium condition in the o-y-direction is identically satisfied. From 
the fact that K, K, lies on the central axis of the bar it follows (Fig. 13) 


1 b  x(D,) b 


Re 
ye ; ag 08 (f+ @,) =F sin (F + ®,). (35) 


+ 


o Copper 
a Aluminium 


Equilibrium along NH FS in the x-direction leads to 


Xds b R c 
aif Ss ; ele —l +2) cos OD, eae +2 ,)sin ®,| 
NH 


Seale +2 @,) 0s ®, nts 1 +2 @,)sin ,| (36) 


| From equ. (35) and (36) R/a and c/a can be found since ®, is known from (34). The restriction 
that c/a = 0 defines the range in wjhich this solution is valid. 
Since the L. H. S. of (36) is always positive it is necessary that 


a eo +2 @,) cos 2,—|5 + 1 +2 ©,)sin ®,| >0 


a y2 |\2 
eo or 
é +1429, 
tan @, < = (37) 
. Ft1+26, 
| or 


DRS Somer 

O07 =. 21-895 (73° 30 Sy = 307). 

| Decreasing y the field remains unchanged, if 2b/a is kept constant, and similar to that for rect- 
_ angular notches. ; 

: 2) Constraint factor. By numerical integration of the moment of the stresses distributed along 
| OANHFS the constraint factor has been calculated and plotted in Fig. 14. For y = 50° and 
| various values of 2 b/a the constraint factor with experimental points are shown in Fig. 15. 

5 
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c) Lower bound. The lower bound is equal to that calculated for rectangular notches (equ. 27) 
wherever y < 7/2 = f, otherwise the side of the stress trapezoid coinsides with the oblique face of 
the notch, 8 = 2/2 —y, and the lower bound of the constraint factor is given by 


Lit Stee =A 6 oe 2 pace 7) (1 = cos y) : (38) 


6. Circular notches. An upper bound for circularly notched bars has been found by Green!. 
A lower bound is shown in Fig. 16: 


t) 1 r /l1—cosp 
a 2 “al cos ). 


Substituting into equ. (21) and (22) we get 
1 iD 2 : 
Dicer =4 lz —* = (1 + sin ) (39) 


a 


and 


(z+ Z)cosb—FeostB—T=0 or PS 15°. (10) 


Upper and lower bounds versus a/(a +r) are plotted in 
Fig. 16. Fig. 17. 


7. V-notched bars. An upper bound solution is due to Green!. The solution is the limiting case 
for trapezoidal notches having 2ba = 0 or for V-notches with circular fillets having a/(a + r) = 0. 
Similarly the lower bound for trapezoidal notches in the limiting case where 2b/a = 0 holds for 
V-notches (Fig. 18). 


Upper Bound (Green) 


2» Aluminium 


/uminium 


$I°' 80° 70° - 66°" 50° 4° 50° = 20° 0 
Fig. 17. Fig. 18. 


8. V-notches with circular fillets. a) Upper bound, Solution I. 1, Slip-line field. The field 
consists of logarithmic spirals in the region Z’OZ.A (Fig. 19), having polar equations 


d = rie”, p=>-?- (41) 


The region AZEF consists of one family of straight lines and their orthogonal trajectories whilst 
the region ZLE, Z’ L’ E’ and AFF’ N are subjected to uniaxial tension. At the newtral point V 


1 A. P. Green, see footnote 3 page 55. 


) 


ot a atl aed 
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there is a stress discontinuity. From equ. 1 it follows 


~~ % + 2 

pend 

_ The unknown quantities to fix the field are y determined from equ. (42) and c/a. From the geo- 
metry of Fig. 19 it follows 


== fo. Ours (42) 


pe ee Ge : 

é a” 2/2. a 2 ; a ee) 
In order to be c/a = 0, 
ee 1 


— oS, . ° 
() 


| Therefore the present solution is valid in that region of the r/a, y plane where both inequalities (42) 
and (44) are satisfied. 


(44) 


Fig. 19. 


2) Constraint factor. The bending moment is equal to twice the moments of the stresses distri- 
buted along OA N with respect to N. By taking these moments the final expression for the con- 
straint factor is 


Lupper = (1 +p) — 52 (e” —1—y) +25 (8? + 2y—4er +3). (45) 


b) Upper bound, Solution II. 1) Slip-line field. When y < (a + 2)/4 = 13.6° a central 
hinge is developed, (Fig. 20). From Henckys theorem applied along BA and ACFS and the fact 
that p, — 0 it follows 


m— 2 AA 


=16.35°, y=y—y- (46) 


Vi 


rf To fix the field the two unknowns must be evaluated. From the geometry of the field, (Fig. 20), 
_ it follows A 
4 Yi Ys 


ba c M4 1+ 2y, Yemen eg herd ee Sets : Ye Vi\ 
3 soe (gt i )+ pein (gt 4 )=3 a\° cos [——— 1]. (47) 


_ Another equation is found from the condition of equilibrium along ACF'S in the x-direction. 
_ (The equilibrium in the y-direction is identically satisfied.) This leads to 


sfee(S + 28) (p+ uals 
1+ 24, : 1+2y, 
ot RE omg) 2) 


(3 Yi + Ve 

be | Saar eee (aE Bee yn Yt 48) 

. = 7e 2 (E+. +) sin ( 5 i: ( 
as 
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From equ. (48) and (49) c/a and R/a are found to fix the field. The coefficient in the equ. (48) is 
always positive. The coefficient of R/a must therefore be positive as well, and this is satisfied for 


ane ea or =p, $ 60°, y 230". 


For y < 30° the slip-line field does not extend to the plane surface of the notch, c/a = 0 and 
solution II for circular notches holds. 


Cipeneeahe ine Riis Da oat I. Nie me Ue 
Fig. 21. 


2) Constraint factor. Taking the moments with respect to the center K, the constraint factor 4 
is evaluated and plotted in Fig. 21. For y = 60° and various values of a/(a + r) the constraint — 


factor together with experimental points is plotted in Fig. 22. 


ee ee a ae 


f 
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c) Lower bound. The lower bound found for circular notches is valid in the present case for 


the same values of a/(a +r) if ySa/2—f. For y = x/2 — B the side of the stress trapezoid 
coinsides with the plane of the notch. In this case: 


p =f-y, Saree =e ar Z | (1 zi cos y) u (49) 


sin y 


a a 


9. Experimental confirmation of the theory. A series of tests was carried out to confirm the 


above theory. Test pieces of high conductivity copper and commertially pure aluminium cold 
_ worked to 30% reduction, were used. Figs. 23 and 24 show the types of stress-strain curves ob- 


- tained. The metal showed a rapid change from the elastic to the plastic condition with a low work 
hardening rate. The yield stress in simple tension was taken as: a) For copper = 22.8 tons 
_ per sq. in., b) for aluminium Y = 8.9 tons per sq. in. The von Mises yield criterion, which is a good 


A 


. 
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approximation for these two metals’, was assumed. Both values of Y agree with the secs: a 
found by plane strain tests, for 30% reduction. The dimensions of the bars used were ‘ 

for copper specimens and 3” x 1?/,”" for aluminium. The notches were sufficiently deep that-a) plane 
strain conditions were insured b) the plastic zone did not spread to the un-notched part of the 


surface. 


The specimens were subjected to four-point loading to give pure bending, using a 50 ton Denison 
s were placed between two pairs of steel rollers and the deflections 


testing machine. The specimen 
h respect 


were measured by dial gauge. Care was taken to secure central positioning of the load wit 
to the axis of the notch, the rollers being well away from the notched region. The moments and 


deflections could be accurately derived from the measurements. 


Tons 
— Stress 


Aluminium 


Strain 
57. 


Fig. 24. 


As soon as the limit moment was reached large plastic deflections occured. The load remained — 
constant but the spring of the testing machine caused instability. To overcome this difficulty the — 
load was increased a little above the desired value to cause rapid rotation of the indicator of the 
gauges and then left to decrease slowly until the rotation stoped. For deflectons well in the plastic 
region, however, continuous readings from both gauges were taken under constant load. By 
plotting the moment-deflection curve and extrapolating back to the elastic line the limit moment My 
was found. Test pieces of length 1 ft were cut from 3 in wide bars of copper and aluminium. Care | 
was taken to secure the accuracy of the dimensions of the notches. A series of V. cutters having | 
semi-angles 30°, 45°, 50°, 55°, 60°, 70° and 75° was used for the V-, Trapezoidal, and V-notches — 
with circular fillets, together with radius cutters R 1/8’, 1/4’, 3/8”, 1/2’ and 1”. For each case © 


Table 1. Experimental results for Rectangular Notches 


Constraint factor 


Dimensions My (tons/in.) 
Symbol | Metal % 
aed ne as ne Theo. Exp. Theo. Exp. toe 
R1 | Cop. 0625 988 1063 3.09 8.79 8.70 1.25 1.238 Ino 
R2 | Al. 125 8163 153 2.975 6.34 6.30 1.217 1.210 —==.0 
R3 | Al. 125 990 227 2.99 2.74 2.71 1.175 1.164 1.0 
R4 | Al. 250 810 309 2.992 2.71 9.65 1.136 1.123 I! 
R5 | Cop. .250 -641 390 3.00 8.96 8.9 1.103 1.096 = .09 
R6 | Al. 375 -736 495 2.985 4,98 4.65 1.067 1.060 — .65 


1G. Lianis and Hugh Ford, J. Mech. Phys. Solids 5 (1956) p. 215. 
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Table 2. Experimental results for Trapezoidal Notches. (Semi-angle y = 50°) 


SS 
Se 


Dimensions My (tons/in.) Constraint Factor 
Symbol] Metal —— 
2b a 26 w 
in. in. Par on Theor. Exp. Theor. Exp. Bye 
DO SS 
Tl Cop. 065 .089 106 | 3.104 8.81 8.80 1.245 1.245 0 
T 2 Al. 125 835 .150 2.995 6.45 6.40 1,205 1.195 — 9 
T3 Al, ALAS) .610 .205 2.980 3.34 3.30 LSet 1.16 —1.0 
T4 Al. .290 .810 .309 2.975 5.64 5.09 1.126 1.11 —1.4 
T5 Cop. 250 -625 400 3.85 8.67 8.65 1.096 1.093 — 3 
T6 Cop. .290 .020 481 | 3.090 5.88 5.83 1.07 1.06 —1.0 
Table 3. Experimental results for Circular Notches 
Dimensions My (tons/in.) Constraint Factor 
Symbol} Metal 
ie i = - ; in Theor. Exp. heor! Exper. ae 
ee ee ee ee 
Cl Cop. .65 61 900 2.905 8.00 7.94 1330 1.320 — .85 
C2 | Cop. 20 -1045 850 3.05 13.0 12.9 1.310 1.30 — .80 
C3 | Al. .250 -7465 . 749 2.930 9.30 5.29 1.263 1.250 —1.0 
C4-} Al. 3125 -970 -645 2.97 2.98 2.93 1.216 1.197 —1.6 
C5 | Al. 500 F023 2000 2.975 2.76 2.73 1.173 1.162 — .95 
C6 | Cop. 900 =| = .426 -460 3.082 4.17 4.15 1.131 1.121 — .90 
Table 4. Experimental results for V-Notches 
Weniensions My (tons/in.) Constraint Factor 
Symbol Metal = 
ae se ee Theor. Exper. Theor. Exper. <i Mie 
| 
Viel: Al. 30° lens 096 3.015 3.80 3.80 1.38 23S 0 
V2 Al. 50° Isa 5 35 3.015 3.02 3.02 1.374 1.374 | 0-7 
V3 Al, 60°.) .584 3.00 3.98 3.57 1.361 1.361 | 0 
V4 Al. G52 a) 000 3.00 3.06 3.63 1.348 1.340 | —.6 
V5 Al. 70° .092 3.00 5.06 3.54 1.320 1.312 —.6 
V6 Al. Toe .088 3.015 3.38 3.36 1.260 1.254 | —.50 
Table 5. Experimental results for V-Notches with Circular Fillets. (y = 60°) 
ee My (tons/in.) Constraint Factor 
Symbol | Metal = = 
oh Ee aE ae Theor. Exper. Theor. Exper. % 
in. in. ot+r in. ly 
cV1] Al. sL2o .817 .871 2.960 6.70 6.70 1.318 1.318 0 
CV 2 Al. 125 485 .195 3.005 2.33 2.31 1.280 1.270, — 8 
CV 3 Al. .250 .629 AAs 3.015 3.81 3.86 1.246 1.241 — A 
CcV4] Al. .2950 sen! -680 2.99 2.67 2.64 1.230 1.218 —1.0 
CV 5 Al. 3125 AT5 .603 3.00 2.08 2.08 1.197 1.197 0 
 CV6| Al. 3125 381 ool 3.015 1.32 1.32 1.173 iL LEGG 0 


six to seven specimens were prepared. The exact notch dimensions were measured by a travelling 
microscope. 


Tables 1 to 5 give experimental details. For tets on Trapezoidal and fillet notches the angle 
of the notch was kept constant and the ratios 2b/a and a/(a + r) were varied, since speciments 
with constant ratios and varying angles would be difficult to prepare to close tolerances. 


: oe hts: it ered pats 
Ss BP ine 
3 
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The experimental results for the constrain factor are given in Figs. 10, 15, 17, 18 and 22. In all 
cases the experimental points are near to the upper bound. The discrepancies are less than 1.5 Yor : 


The results from these tests show that the upper bound is the complete solution, and the fact 
that the experimental points lie slightly below the upper bound, is due to a slight overestimation | 
of the yield stress. 

10. Conclusion. It is believed that all practical cases of notch shapes are covered by the above 
analysis. The experiments are in favour of the upper bound solution, though, for practical cases, — 
the mean value derived from the limit analysis differs very little from that of the bounds, and it © 


can be used with confidence. 4 


Further tests, however, are needed to confirm the deformation pattern. It is also believed — 
that this analysis will furnish useful informations on the study on the mechanism of fracture. 


(Eingegangen am 5. Mai 1959.) 


Anschrift des Verfassers: Professor George Lianis, Ph. D., D. I. C., Dipl.-Ing., 
School of Aeronautical Engineering, Purdue University, Lafayette, Indiana, U.S.A. 
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Stufenlos verstellbare mechanische Getriebe 


Von Dipl.-Ing. FrrppR. W. Smmonis VDI, Reg.-Baurat a. D., Beratender Ingenieur VBI, 
Vereidigter Sachverstandiger fiir Werkzeugmaschinen und Werkzeuge, Berlin 


Zweite, vollig neubearbeitete und erweiterte Auflage 
des in erster Auflage erschienenen Werkstattbuches Heft 96 
Mit 252 Abbildungen. VIII, 190 Seiten Gr.-8°. 1959. Ganzleinen DM 29,40 


INHALTSUBERSICHT 


Einleitung - Reibradgetriebe mit Kegelscheiben und Zylindern - Reibradgetriebe mit umlaufendem Reib- 
ring - Reibgetriebe unter Benutzung von Kugeln als Ubertragungsmittel (Cloboid-Getriebe) - Planeten- 
getriebe mit frei umlaufenden Kugeln oder Rollen - Umhiillungs- oder Umschlingungsgetriebe mit Riemen- 
ubertragung - Umbiillungs- oder Umschlingungsgetriebe mit Ketteniibertragung - Theorie der Umschlin- 
gungsgetriebe mit keilformigen Umlaufflachen - Schaltwerksgetriebe - Steigerungsméglichkeiten fir iiber- 
tragbare Leistung und erzielbaren Verstellbereich bei stufenlosen Getrieben (Differentialgetriebe) - Selbst- 
tatige Steuerungen von stufenlosen Getrieben zur Erfiillung eines geforderten bestimmten Arbeitsablaufes - 
Sondergetriebe - Sachverzeichnis. 


Optimale Stufenridergetriebe fiir Werkzeugmaschinen 


Errechnung und Raderanordnung 


Von Dr.-Ing. E.h. ERWIN STEPHAN 
Mit 139 Abbildungen. VI, 102 Seiten Gr.-8°. 1958. DM 21,—; Ganzleinen DM 25,— 


AUS DEN BESPRECHUNGEN 


»,Das Buch soll bei einem vorgeschriebenen Drehzahlbereich und fiir eine vorgegebene Stufenzahl Wege 
zeigen, die das Auffinden optimaler Stufenradergetriebe erleichtern. Fiir neuzeitliche Anspriiche ist es 
erwiinscht, einem solchen Getriebe eine geometrische Drehzahlstufung zu geben, Durch die Zweiteilung in 
die Abschnitte ,Entwurf und Errechnung‘ und ,Bauliche MaSnahmen‘ tragt der Verfasser fiir den Entwurf 
derartiger Geiriebe Rechnung. Wegen der Méglichkeit, Aufwand und Platzbedarf von Stufenradergetrieben 
durch Anwenden gebundener Getriebe herabzusetzen, geht er auf die einfache und die mehrfache Bindung 
besonders ein. Auf die Schwierigkeiten, da& bei doppelt gebundenen Getrieben die Normdrehzahlen nur 
schwer und bei dreifach gebundenen nicht einzuhalten sind, weist der Autor besonders hin. Von grofem 
Wert ist die Beschreibung der Einschaltméglichkeiten der Rader wahrend des Laufes mit Hilfe elektrisch 
oder hydraulisch betatigter Reibungskupplungen bzw. durch Kupplungen mit Synchronisiereinrichtungen. 
Das kleine iibersichtliche Werk kann dem einschlagigen Konstrukteur manche Anregung vermitteln.‘‘ 
VDI-Zeitschrift 
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Kinematisch-getriebedynamisches Praktikum 


Lehr- und Ubungsbuch zur graphodynamischen Analyse ebener Getriebe 
fiir den Konstrukteur, die Vorlesung und das Selbststudium 


Von Dr. phil. habil. RUDOLF BEYER, apl. Professor fiir Getriebelehre und Kinematik an der Tech- 
nischen Hochschule Miinchen, Oberstudienrat a. D. des Oskar-v.-Miller-Polytechnikums, Akademie fiir 
angewandte Technik, Miinchen 


Mit 125 Abbildungen. VIII, 170 Seiten Gr.-8°. 1960. Ganzleinen DM 29,40 


INHALTSUBERSICHT 


Grundlagen. Kraftreduktion in ebenen Getrieben. Kinetische Energie des Cetriebes. Das d’Alembertsche 
Prinzip. Die dynamische Grundgleichung fiir die Kraft- und Massenwirkungen in zwanglaufigen Cetrieben. 
Beispiele zur I. Wittenbauerschen Grundaufgabe. Zusammenfassendes Beispiel. Kraft- und Massenwir- 
kungen an einem siebengelenkigen Koppel-Rastgetriebe. Korrektur der Stabkrafte. Das Federhoferschen 
A-Verfahren. Dynamik der Getriebe mit doppeltem oder mehrfachem Antrieb. Erganzende Grundlagen 
zur Dynamik des komplan bewegten Getriebegliedes. Bewegung des Schwerpunktes der Gliederanordnung 
eines ebenen Getriebes. Das Wittenbauersche Massen-Wucht-Diagramm. Reibung in Getrieben. Ubungs- 
beispiele mit Lésungen. Aufgaben. Schrifttum. Namen- und Sachverzeichnis. 


Vom gleichen Verfasser erschienen: 


Kinematisch-getriebeanalytisches Praktikum 


Hand- und Ubungsbuch zur Analyse ebener Getriebe 
fiir den Konstrukteur, die Vorlesung und das Selbststudium 


Mit 162 Abbildungen. WII, 172 Seiten Gr.-8°. 1958. Ganzleinen DM 28,50 


Kinematische Getriebesynthese 


Grundlagen einer quantitativen Getriebelehre ebener Getriebe 
fiir den Konstrukteur, fiir die Vorlesung und das Selbststudium 


Mit 258 Abbildungen. VII, 217 Seiten Gr.-8°. 1953. Ganzleinen DM 36,— 
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